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ABSTRACT 


A  generally  applicable  method  of  obtaining  kinetic  parameters  from 
temperature-programmed  thermogravimetry  is  presented.  Fetors  influencing 
the  selection  of  a  particular  method  for  the  numerous  treatments  reported  in 
the  literature  are  discussed  in  detail.  The  method  of  Friedman  involving  the 
use  of  several  thermograms  at  different  heating  rates  and  determining  Arhennius 
parameters  at  each  pe:  ant  conversion  was  chosen.  The  experimental  procedure 
and  a  method  of  handling  thermogravimetric  analysis  (TGA)  data  and  calculations 
by  computer  are  fully  described.  The  application  of  the  treatment  to  some 
specific  polymer  degradation  systems  is  reported  in  order  to  illustrate  the 
scope  of  the  method  and  its  potential  usefulness  in  obtaining  information  con¬ 
cerning  complex  degradation  mechanisms.  Poly(tetrafluoroethylene),  an  ali¬ 
phatic,  and  an  aromatic  polyamide  were  the  polymers  selected  for  this  study. 
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SECTION  I 


INTRODUCTION 


A  knowledge  and  understanding  of  polymer  degradation  processes  is  nec¬ 
essary  to  improve  the  performance  of  polymers  in  high  temperature  applications 
and  to  direct  research  towards  more  thermally  stable  systems.  One  of  the 
most  important  methods  of  studying  polymer  degradation  is  to  examine  the 
kinetics  governing  the  breakdown  reactions  from  which  information  about  the 
thermal  stability  and  the  mechanism  of  degradation  can  be  obtained. 

In  degradation  kinetics,  the  rate  of  change  of  some  property  must  be  mea¬ 
sured  as  a  function  of  time,  temperature,  or  conversion;  any  of  several  proper¬ 
ties  can  be  utilized  if  the  property  can  be  correlated  with  the  degradation 
process,  e.g.,  mass,  molecular  weight,  a  chemical  group,  or  production  of 
volatile  products.  A  common  method  has  been  to  follow  weight  changes.  A 
thermobalance  can  be  used  to  continuously  record  the  weight  changes.  In  the 
past,  isothermal  studies  have  been  most  common,  in  which  the  sample  is  main¬ 
tained  at  a  constant  temperature  while  weight-time  measurements  are  made.  In 
recent  years,  temperature-programmed  methods  have  been  increasingly  used. 
In  this  method,  temperature  is  continuously  raised,  usually  linearly  with  time, 
and  a  thermogram  of  weight  versus  temperature  obtained.  In  theory,  the  pro¬ 
grammed  method  should  have  certain  advantages  over  the  isothermal  method. 
An  important  criticism  of  isothermal  methods  is  that  the  temperature  of  the 
sample  cannot  be  raised  instantaneously  to  the  desired  temperature  so  that 
some  weight  may  be  lost  before  the  degradation  temperature  is  reached.  This 
weight  loss  may  give  rise  to  an  apparent  maximum  in  the  rate  curve  which  masks 
the  true  initial  features.  Initial  rate  characteristics  are  very  important  since 
end-group  and  impurity- induced  reactions  may  show  up.  Temperature- 
programmed  methods  should  overcome  this  problem.  It  was  further  visualized 
that  since  weight-temperature  dependencies  were  contained  in  a  single  chart 
the  thermogram  from  a  programmed  experiment  would  yield  the  equivalent 
information  of  a  large  family  of  isothermal  experiments.  Thus,  sample  uni¬ 
formity  errors  would  be  avoided  and  the  procedure  would  be  less  time  consuming. 
The  thermograms  are  also  useful  for  qualitative  comparisons  of  thermal 
stabilities. 
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A  =  pre-exponential  factor 
Ea  =  activation  energy 

R  =  general  gas  constant 
T  =  absolute  temperature 

The  normal  kinetic  approach  is  to  determine  rate  constants  for  a  given  reaction 
at  various  temperatures  using  Equation  1  or  an  integrated  form  of  it  and  lence 
to  calculate  the  parameters  A  and  E  from  Equation  2.  The  activation  energy 

a 

can  often  be  correlated  with  the  breaking  of  specific  chemical  bonds  and  gives 
important  information  concerning  the  mechanism  of  the  reaction  being  studied. 

Polymer  degradation  kinetics  are  normally  studied  in  the  solid  or  melt  phase 
and,  since  a  chemical  reaction  is  occurring,  it  is  assumed  that  a  kinetic  treat¬ 
ment  and  rate  expression  is  applicable.  Owing  to  the  complexity  and  variety  of 
polymer  decomposition  schemes,  it  is  found,  however,  that  the  concentration 
of  polymer  molecules  is  not  equivalent  to  the  concentration  of  reactant  in 
normal  reactions  so  that  the  term  “concentration”  must  be  used  with  extreme 
caution  in  discussing  polymer  degradation.  Two  well  established  types  of  deg¬ 
radation  illustrate  this  point  (Reference  1).  In  the  first,  degradation  is  initiated 
at  a  chain  end  and  proceeds  by  unzipping  through  the  entire  polymer  chain.  The 
sample  loses  one  polymer  molecule  and  its  corresponding  weight  and  volume, 
but  the  concentration  of  polymer  molecules  remains  unchanged.  In  the  second 
type,  the  elimination  of  side-groups  in  the  chain  alters  the  chemical  structure 
of  the  polymer  but  does  not  necessarily  influence  the  number  or  concentration 
of  polymer  molecules.  In  polymers,  discussion  in  terms  of  reactive  sites  which 
could  be  chain  ends,  particular  bonds  or  groups  in  the  structure,  the  links 
Joining  monomer  units,  etc.,  is  often  more  fruitful  for  elucidation  of  mechanisms. 

ft  is  found  tn  practice  that  most  degradation  reactions  result  in  loss  of 
weight  owing  to  the  formation  of  smaller,  volatile  species  ^  rearrangements  such 
as  the  coloration  of  poly(acrylonitrlle)  are  notable  exceptions  J.  Therefore,  it 
has  been  found  convenient  to  study  kinetics  in  terms  of  weight  loss;  to  this 
extent,  the  approach  is  empirical  since  the  only  proof  of  validity  is  that  exper¬ 
imental  results  appear  to  fit  the  adopted  kinetic  expressions.  In  order  to 
preserve  a  close  resemblance  to  normal  kinetic  procedure,  the  weight  terms 
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should  be  expressed  in  fractional  form.  Further,  since  a  reactant  should  have 
aero  concentration  on  completion  of  reaction,  the  weight  term  should  be  cor¬ 
rected  for  any  residue  weight  remaining  after  degradation.  A  general  rate 
expression  based  on  weight  terms  can  now  be  postulated 


I  dW  /  W-W<  \ 

%  ~  Wf  df  "  '  wo-wf  '  = 

W  =  instantaneous  weight  of  sample 
WQ  =  initial  weight  of  sample 

=  final  or  residue  weight 


k  F(W) 


W  -Wf 

rr. — ttj—  =  fractional  weight  remaining 
Wo"wf 


W  -Wf  WQ-W 

(1  “  w"w  "  =  w°-w  ■  =  fractional  weight  loss) 
o  f  o  f 


(3) 


Equation  3  is  equivalent  to  that  used  by  Doyle  (Reference  2)  except  that  weight 
terms  are  retained  here  in  preference  to  Doyle’s  active  weight  fraction,  h.  It 
must  be  emphasized  that  the  functionality  F(W)  need  not  be  simple  and  should 
not  be  assumed  in  advance.  The  presumption  of  “order”  type  functionality  should 
be  avoided  unless  it  is  warranted  by  experimental  evidence.  Two  particular  forms 
of  F(W),  which  have  beeu  observed  in  isothermal  studies,  merit  a  brief  dis¬ 
cussion  at  this  point: 


a.  F(W)  =  yj—  (for  simplicity,  assume  Wf  =  O) 
o 


This  is  a  common  cr "e  which  results  in  apparent  first  order  kinetics  by 
analogy  to  classical  rate  laws.  The  rate  of  weight  loss  is  proportional  to  the 
residual  weight  of  polymer.  Polymers  which  degrade  predominantly  by  unzipping 
show  kinetics  of  this  type,  the  weight  loss  being  a  direct  measure  of  the  amount 
of  unzipping  and  degree  of  degradation.  This,  of  course,  is  a  simplified  view  of 
the  overall  kinetics,  and  many  other  interesting  kinetic  dependencies  can  be 
observed  (References  1  and  3). 


b 


F(w)  =  (  L  -  i ) 


L  N  -  L  )  L  -  I  ) 

a )  ■  - 

N 


(4) 
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L  -  smallest  chain  length  which  does  not  evaporate  under  degradative 
conditions 

N  =  initial  chain  length  of  polymer 
-kt 

a  =  1-e  (temperature  constant)  =  degree  of  degradation  (bond  scission) 

This  complex  form  of  F(W)  has  been  developed  by  Simha  et  al  as  the  solution  to 
the  random  type  of  degradation  mechanism  (References  4  and  5).  This  mechanism 
consists  of  chain  rupture  taking  place  randomly;  the  act  of  scission  need  not 
result  in  weight  loss  which  is  possible  only  when  scission  occurs  near  a  chain 
end  to  give  a  chain  fragment  of  length  <L.  The  rate  of  weight  loss,  corresponding 
to  a  spectrum  of  products  of  chain  length  from  one  to  L,  is  proportional  to  the 
number  of  chain  ends  in  the  system,  which  increases  at  first  owing  to  the  random 
scission  of  the  chains „  Eventually,  the  average  chain  length  becomes  small 
enough  that  further  increase  in  the  number  of  ends  is  balanced  by  loss  through 
evaporation.  The  rate  of  weight  loss  passes  through  a  maximum  and  subsequently 
falls ,  it  has  been  shown  theoretically  that  the  maximum  should  occur  at  ap¬ 
proximately  26%  conversion.  Polymethylene  and  polypnethyl  acrylate)  are 
examples  of  vinyl  polymers  exhibiting  the  above  features.  Many  condensation 
type  polymers  such  as  cellulosics,  polyamides,  and  polyesters  also  show  random 
characteristics  (Reference  3). 

As  has  been  discissed,  several  functionalities  of  weight  are  readily  ob¬ 
served  in  isothermal  weight-loss  studies.  In  contrast,  primary  data  from 
programmed  TGA  studies  is  much  less  revealing  nd  differences  caused  by 
changes  in  weight  functionality  may  be  very  subtle.  The  weight-temperature 
thermogram  and  the  derivative  rate-temperature  curve  (Figure  1)  are  somewhat 
general  in  character,  apart  from  the  obvious  difference  between  a  simple  and 
multistep  process;  differences  in  shape  are  relative  and  could  be  divined  only 
by  reference  to  other  “known”  curves.  That  both  activation  energy  and  mech¬ 
anism  affect  the  shape  of  the  curve  adds  to  the  difficulties  (Reference  6).  As  is 
shown  in  this  reference,  rate-conversion  curves  are  of  greater  diagnostic 
importance,  although  in  practice,  experimental  fluctuations  from  a  smooth  curve 
might  easily  mask  the  diagnostic  features,  namely  the  conversion  for  the  max¬ 
imum  rate  and  the  initial  slope  of  the  curve.  The  ability  of  kinetic  analysis  to 
afford  not  only  the  Arrhenius  parameters  but  also  supplementary  information 
about  the  mechanism  is  investigated  in  the  following  paragraphs. 
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2.  KINETIC  ANALYSIS  OF  PROGRAMMED  WEIGHT-LOSS  DATA 

Analytical  methods  can  be  broadly  divided  into  those  using  differential  and 
Integral  treatments.  Differential  methods  are  based  on  Equation  3.  In  pro¬ 
grammed  TGA,  k  In  this  equation  becomes  a  variable  since  temperature  and 
therefore  the  rate  conotant  are  continuously  changing.  The  expression  can  be 
rearranged  as  shown  for  a  linear  heating  rate,  B, 


B 

I _  dW^ 

VV  dT 


dT 

dt 

A  -E0/RT 

T  '  ® 


F(W) 


(5) 


The  integral  approach  was  developed  to  use  TGA  weight-loss  data  directly  and 
to  avoid  calculation  of  rates.  The  equation  of  the  thermogram  (References  2 
and  7)  is  obtained  by  integrating  Equation  5 


dW 


(W0-Wf)F{W) 


,  A  r' 

-  F  =  "F  •  / 
T0 


exp 


The  result  is  given  as 


F'(W)  « 


HE  q 
BR 


RT  / 


(6) 


(7) 


in  which  p(E  /RT)  is  a  complex  integral,  hi  practice,  Integral  methods  use 
some  convenient  approximation  based  on  Equation  7. 


A  detailed  critique  of  all  the  previously  reported  kinetic  treaments  is 
beyond  the  soope  of  this  report.  Several  reviews  have  been  published  (References 
2,  8,  aud  9)  and  the  reoent  paper  by  Flynn  and  Wall  (Reference  6)  offers  an 
excellent  critical  survey  of  many  of  the  methods.  It  is  fitting,  however,  to 
examine  some  of  the  approaches  and  evaluate  their  usefulness  for  obtaining 
meaningful  information  about  the  degradation.  Some  representative  approaches 
are  listed  in  Table  L  It  is  found  that  almost  all  the  methods  are  subject  to  one 
or  more  of  the  following  criticisms: 


a.  The  equation  fits  a  limited  portion  of  the  curve  so  that  parameters  for 
the  whole  degradation  range  ca.mot  be  obtained. 
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b.  Degradation  kinetics  axe  assumed  to  be  of  the  “order”  type,  and  per¬ 
haps  even  first  order  kinetics  are  specified.  Since,  as  has  been  previously 
discussed,  this  assumption  is  not  necessarily  valid,  if  is  impossible  to  conclude 
whether  the  resulting  parameters  are  meaningful.  Often  a  random  type  deg¬ 
radation  and  probably  other  types  will  give  a  deceptively  good  fit  to  the  first 
order  approximation  over  most  of  the  conversion  range.  Errors  in  assumed 
order,  however,  give  vast  errors  in  the  other  kinetic  parameters;  thus,  in 
Reference  6,  a  random  curve  corresponding  to  an  activation  energy  of  60 
kcal-mole-1  was  analyzed  by  a  first  order  treatment  to  have  an  activation  energy 
ranging  from  66  to  175  kcal-mole-1  depending  on  conversion. 

c.  The  time  saving  advantage  of  single  point  and  single  curve  methods  is 
obvioua^ut  it  cannot  be  overemphasized  that  such  analyses,  when  applied  to  a 
new  s$*em  of  unknown  kinetics,  may  be  absolutely  worthless  with  regard  to 
significant  kinetic  parameters. 

d.  Care  must  be  taken  when  deriving  approximate  solutions  since  a  sequence 
of  steps  Involving  successive  approximations  may  give  a  final  working  equation 
which  is  not  particularly  applicable.  Assumptions  and  approximations  should 
be  kept  to  a  minimum. 

On  this  basis,  methods  involving  a  single  thermogram  can  be  discarded  and 
the  potential  of  methods  using  several  thermograms  examined.  The  method  of 
Reich  using  two  heating  rates  (Table  I)  is  ^satisfactory  since  no  information 
about  uie  rate  law  can  be  obtained.  Chatter jee’s  method  involving  different 
sample  weights  must  also  be  discarded  since,  in  addition  to  the  assumption  of 
order,  the  reasoning  is  erroneous  (Reference  37).  When  weight  terms  are 
expressed  in  the  proper  fractional  units,  his  treatment  collapses  completely.  In 
the  equation  for  order  (Table  I),  it  can  be  visualized  that  “n”  would  always  be 
unity,  irrespective  of  the  mechanism  of  degradation,  excepting  the  possible 
influence  of  experimental  irreproducibllity  and  diffusion  controlled  weight  loss. 

Therefore,  the  choice  of  methods  is  reduced  to  two,  one  integral  and  the 
other  differential;  both  require  a  series  of  thermograms  to  be  obtained  using 
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different  heating  rates  and  appear  to  be  similar  in  potential  and  general  appli¬ 
cability.  This  also  represents  the  conclusion  of  the  Flynn  and  Wall  survey.  In 
this  report,  selectivity  was  carried  one  stage  further  according  to  the  following 
arguments: 

a.  hi  the  past,  integral  methods  have  had  the  advantage  that  tedious  manual 
derivations  of  rates  of  weight  loss  could  be  avoided.  Since  many  thermo¬ 
balances  now  incorporate  automatic  differentiation  which  allows  simultaneous 
recording  of  TGA  and  differential  thermogravimetry  pTG)  curves  and  since, 
lacking  this  equipment,  rates  can  be  computed  from  weight-loss  data  (see 
Section  IQ),  differential  methods  need  no  longer  be  avoided. 

b.  The  integral  treatment  involves  more  approximations  than  the  differential 
method.  Thus,  the  complex  integral  is  normally  approximated  to  a  convenient 
working  equation. 

c.  The  integral  treatments  assume  a  linear  heating  rate  throughout  the 
degradation.  This  assumption  is  not  strictly  valid  since,  for  example,  volatil¬ 
ization  of  gaseous  products  produces  a  cooling  effect  on  the  polymer  and  a 
lowering  of  the  heating  rate. 

d.  In  order  to  determine  the  rate  law  dependence,  F(W),  in  the  integral 
method,  experimental  curves  must  be  fitted  to  theoretical  curves  for  known 
kinetic  schemes.  Experimental  errors  and  kinetic  irregularities  (the  latter  are 
often  observed  during  the  initial  and  final  stages  of  weight  loss)  would  create 
difficulties  in  curve  fitting.  A  more  serious  disadvantage  is  that  the  true 
kinetic  scheme  might  be  different  from  the  available  theoretical  ones,  e.g.,  the 
isothermal  rate  curve  maximum  at  40%  con  vers  ion  for  polystyrene  (Reference  22). 
hi  the  study  of  the  complex  condensation  polymers  currently  being  used  for  high 
temperature  applications,  new  types  of  rate  dependencies  may  well  appear  and 
would  necessitate  laborious  trial  and  error  curve  fitting  in  the  integral  method. 
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3.  “FRIEDMAN”  METHOD  OF  ANALYSIS 

The  foregoing  arguments  resulted  in  the  selection  of  the  method  first  used 
by  Friedman  (Reference  15).  The  basic  equation  used  is  shown  in  Table  I  and  is 
the  logarithmic  expression  of  Equation  3 

log  l  (W  -  wTT  '  TT  1  =  loflA  +  log  F<W)"  2.303RT  (8) 

0*0  D 

Equation  8  denotes  that  both  the  rate  and  temperature  corresponding  to  any 
specific  value  of  F(W)  are  dependent  on  the  heating  rate  employed.  If  it  is 
assumed,  for  the  moment,  that  the  chemistry  of  the  orocess  is  independent  of 

Wo"W 

temperature,  then  F(W)  will  be  constant  at  any  weight  loss,  — — tjt— »  regardless 

WQ-wf 

of  heating  rate.  Thus,  by  measuring  or  calculating  the  rates  and  temperatures 
corresponding  to  the  same  fractional  weight  loss  in  a  series  of  thermograms 
obtained  at  different  heating  rates  and  by  plotting  log(rate)  against  reciprocal 
temperature,  the  activation  energy  can  be  calculated  for  any  conversion.  An 
examination  of  the  changes  in  calculated  activation  energies  over  the  whole 
conversion  range  reveals  whether  the  kinetics  vary  with  conversion  although  it 
may  be  difficult  to  resolve  experimental  errors  and  real  changes  in  activation 
energy.  If  the  activation  energy  does  not  vary  much  over  the  entire  range,  an 
average  value,  E  ,  is  calculated.  Hence,  by  rearranging  Equation  8 

loqAF(W)  =  Ioq  A  +  log  F(W)  =  log  (rate)  -►  —  —  (9) 

Values  of  log  A  F(\v)  are  calculated  at  each  conversion  for  the  various  heating 
rates.  Theoretically,  there  should  be  no  heating  rate  variation,  A  and  F(W)  being 
considered  to  be  independent  of  temperature  in  the  simple  case,  but  experimental 
errors  usually  result  in  a  small  spread  of  log  A  F(W)  values  so  that  an  average 
value  is  again  calculated. 

When  this  process  has  been  repeated  for  all  conversions,  a  plot  of  log  A  F(W) 
w  W  “ 

against  log  (or  log  yfZw"  for  degradations  producing  a  residue)  reveals 
o  o  f 

the  rate  law.  In  the  programmed  TGA  method,  this  plot  removes  the  effect  of 
heating  rate  and  ohanglng  temperature;  for  diagnostic  purposes,  it  is  analogous 
to  the  rate-conversion  curve  in  isothermal  treatments.  In  many  cases,  therefore, 
the  appropriate  weight  functionality  will  be  elucidated  (Table  II). 
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TABLE  II 

EFFECT  OF  F(W)  ON  LOG  A  F(W)  PLOT 
F(W)  LOG  A  F(W)  PLOT 


W  W  ~W 

:  w -w, 

o  of 

/  W  -W  \n 

\  w -w,  I 

O  I 

“random”  function 


straight  line,  slope  =  1 

straight  line,  slope  =  n 

curve,  maximum  at  log(0.74) 
(26%  conversion) 


Other  types  of  kinetics,  such  as  those  giving  a  rate  maximum  at  conversions 
other  than  26%,  would  be  just  as  easily  resolved  by  this  plot.  Kinetic  irregu¬ 
larities  are  revealed;  for  example,  early  weight  loss  caused  by  lower  activation 
energy  processes  normally  shows  as  a  steeply  falling  portion  in  the  initial 
stages  of  the  log  A  F(W)  curve.  In  fact,  any  true  change  in  the  kinetics  will 
produce  a  change  in  slope  of  the  plot  based  on  the  assumption  of  a  single 
activation  energy  process.  The  extreme  sensitivity  of  log  A  F(W)  to  changes 
in  E  explains  why  an  average  E„  i3  used  instead  of  individual  values  at  each 
conversion.  If  the  latter  were  used,  the  experimental  fluctuations  of  E  would 
outweigh  the  effect  of  rate  on  log  A  F(W);  the  resulting  plot  would  be  very 
scattered  and  no  Information  about  the  weight  functionality  would  be  obtained. 
Other  potential  sources  of  error  in  determining  the  form  of  F(W)  are  wrongly 
assigned  conversion  ranges  and  use  of  wrong  conversion  units.  These  may 
seem  obvious  but  the  former  can  be  easily  done  in  some  complex  thermograms 
and  the  latter  follows  from  perusal  of  the  variety  of  expressions  used  in  the 
past.  For  example,  in  the  original  description  of  this  method,  Friedman  used  a 


W-Wf 

different  expression  of  concentration^  ) ,  from  that  proposed  in  this 

Wo 

report.  The  result  was  a  very  high  apparent  order  of  reaction. 


The  pre-exponential  factor  can  be  easily  calculated  by  subtracting  the 
F(W)  value  from  the  log  A  F(W)  term.  The  determination  of  A  could  be  ac¬ 
complished  graphically  by  replotting  log  A  F(W)  against  log  F(W)  and  ex¬ 
trapolating  the  resulting  straight  line  to  log  F(W)  =  O.  The  ability  of  the  method 
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to  completely  analyze  more  complex  degradation  systems  with  some  typical 
complicating  features  is  discussed  in  the  following  paragraphs. 

a.  Activation  energy  changes  with  conversion 

Certain  degradation  mechanisms  involve  real  activation  energy  changes 
with  conversion  (References  3  and  5).  The  reason  may  be  a  dependence  of  E 

& 

on  the  molecular  weight  which  itself  varies  with  conversion,  a  change  in  mech¬ 
anism  as  in  poly(methyl  methacrylate),  or  a  change  in  structure  of  the  polymer. 
Suppose  a  smooth  increase  in  E  is  observed  in  the  E  -conversion  plot  obtained 
by  the  preceding  analysis.  The  assumption  of  an  average  E  to  calculate  log 

a 

A  F(W)  would  invalidate  the  rate  law  determination.  Features  such  as  rate 
maxima  would  still  be  discerned  but  the  overall  curve  would  be  skewed  relative 
to  the  '‘theoretical”  curve.  Then,  it  is  conceivable  that  a  smooth  E  profile 

fl. 

might  be  used  instead  of  an  average  E  to  give  mors  meaningful  information 

a 

about  F(W). 


b.  Random  Degradations 

The  complex  functionality  of  weight  for  this  mechanism  (Equation  4)  casts 
doubt  on  the  validity  of  the  assumption  in  the  kinetic  analysis  that  F(W)  is 
independent  of  temperature.  The  source  of  this  contention  is  that  a,  the  degree 
of  degradation  inherent  in  the  function,  contains  a  temperature  term.  It  can  be 
shown,  however,  that  a  is  independent  of  heating  rate  and  depends  o.Jy  on  the 
conversion.  Under  dynamic  conditions 

t 

-  J  kdt 

a  =  I  -  e  o 

a  =  I  -  e  B  J  rtB  (10) 

'o 

That  a  is  independent  of  B  follows  from  integral  treatments,  particularly  that 
in  Reference  21.  The  method  is  still  valid  therefore  for  the  random  case.  The 
derivation  of  a  pre-exponential  factor  in  the  random  case  may  be  difficult  since 
the  various  parameters  needed  to  calculate  F(W)maynot  be  available.  Equation  4 
is  an  approximate  solution  dependent  on  certain  boundary  conditions  and 
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may  not  be  applicable  to  some  real  cases  although  the  overall  random  curve 
is  still  observed.  The  normal  method  of  obtaining  A  from 

/  dC  \ 

( — —  =  L/e  (Reference  5)  (II) 

man 

and  in  this  treatment. 

flogAF(W)]  =  log  A  -i  log  (-^t)  (12) 

1  Jmas  '  ' 

must  be  used  cautiously. 

c.  Complex  Mechanisms 

The  treatment  of  two  of  the  more  straightforward  complex  cases  was  dis¬ 
cussed  in  Reference  6.  The  first  case  involved  competitive  reactions  in  which 
the  rate  curve  and  thermogram  appeared  similar  to  that  for  a  simple  reaction 
except  for  irregular  trends  in  the  maximum  rate.  The  ability  of  the  method  to 
resolve  the  two  reactions  depends  on  how  different  the  individual  parameters 
are.  If  they  have  similar  orders  and  activation  energies,  it  is  doubtful  whether 
any  resolution  could  be  achieved.  The  second  case  consisted  of  two  independent 
reactions,  each  of  which  could  be  observed  in  thermograms  obtained  at  low 
heating  rates.  The  corresponding  activation  energies  were  obtained  at  low  and 
high  conversions.  Only  the  two  methods  involving  several  heating  rates  showed 
any  success  in  resolving  theee  cases.  Inreal polymer  degradations,  the  following 
complex  cases  have  been  observed  in  previous  and  current  work: 

a.  The  thermogram  consists  of  several  consecutive  steps  with  distinct 
plateaus  between  the  decompositions.  This  case  Is  easily  dealt  with  by  treating 
each  step  individually  as  a  simple  case. 

b.  The  thermogram  exhibits  overlapping  reactions  and  the  DTG  curve  has 
several  maxima.  This  could  be  visualized  as  Case  a.  in  which  the  second  step 
commences  before  the  first  reaction  is  complete.  This  is  not  strictly  analogous 
to  the  “Independent  reaction”  discussed  In  “Complex  Mechanisms”  since  the 
amount  of  each  reaction  may  depend  on  the  heating  rate.  In  such  a  case,  a  proper 
analysis  may  be  thwarted  although  relevunt  information  would  still  be  obtained 
for  the  low  temperature  reaction  by  studying  the  initial  portion  of  the  weight  loss. 
An  example  of  this  type  Is  presented  In  Section  VL 
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c.  A  complex  curve  somewhat  similar  to  b.  has  been  observed  for  some 
aromatic  polyesters  (Reference  23).  For  these  polymers,  the  major  weight-loss 
reaction  changes  smoothly  into  a  slow  char-forming  reaction  which  gives  the 
rate  curve  a  long,  high  temperature  tail.  Once  again,  it  is  predicted  that  only 
limited  information  will  be  obtained  and  further  discussion  must  await  detailed 
examination  of  actual  examples. 

4.  SIGNIFICANCE  OF  KINETIC  PARAMETERS 

The  influence  of  experimental  variables  in  programmed  thermogravimetry 
and  their  effect  on  the  resulting  kinetic  parameters  has  been  the  subject  of 
several  reviews  (References  2,  8,  19,  24,  and  25)  and  has  convinced  some 
authors  that  the  parameters  are  purely  empirical.  However,  employing  careful 
techniques  and  strict  standardization,  many  of  the  sources  of  error  such  as 
weighing  errors,  diffusion  effects,  and  differences  dependent  on  the  physical 
form  of  the  sample  can  be  minimized  or  eliminated.  The  sources  of  error,  which 
could  be  considered  appropriate  to  programmed  methods  and  not  to  isothermal 
techniques,  are  heating  rate  and  temperature  errors.  The  latter  can  be  removed 
by  good  experimental  procedure  and,  as  has  been  shown  previously,  the  former 
effect  is  removed  in  the  ultimate  analysis  by  the  method  chosen.  It  is  concluded, 
therefore,  that  the  programmed  method  used  should  give  information  equivalent 
to  that  obtained  isothermally  and  it  is  contended  that  this  information,  especially 
the  activation  energy,  can  be  meaningful  with  respect  to  stability  and  mechanism. 
Thus,  as  in  Reference  26,  overall  activation  energies  have  been  related  suc¬ 
cessfully  to  the  energies  of  individual  steps  comprising  the  reaction.  A  better 
proof  must  be  to  compare  the  experimental  activation  energy  to  that  observed  in 
conventional  kinetic  analysis.  This  should  be  feasible  when  a  polymer  and  its 
model  degrade  by  exactly  similar  mechanisms.  The  field  of  condensation 
polymers  may  contain  examples  satisfying  this  condition. 
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SECTION  m 

EXPERIMENTAL  PROCEDURE 

1.  POLYMER  SAMPLES 

Descriptions  of  the  methods  employed  for  the  preparation  and  purification 
of  samples  are  given  in  the  reports  which  describe  in  detail  the  results  for 
those  polymers, 

2.  APPARATUS 

The  thermobalances  used  were  the  Ainsworth  Models  AVand  RV  which  gave 
full  scale  recorder  deflections  of  100  mg  and  10  mg,  respectively.  In  most 
cases,  several  deflections  of  the  recorder  pen  were  necessary  to  follow  the 
complete  weight  loss  of  samples. 

A  sectional  diagram  of  the  thermobalance,  degradation  tube,  and  furnace 
is  given  in  Reference  27. 

The  temperature  programmer  used  was  the  West  Gardsman  Model  JGB 
Program  Controller  which  operated  a  proportioning  power  supply.  The  temper¬ 
ature  set  point  was  driven  by  a  cam  cut  to  give  close  to  a  linear  increase  in 
temperature  with  time.  Variation  in  pi  "am  rate  was  effected  by  changing 
gears  in  the  motor  to  cam  gear  train.  Sorrv  uf  the  approximate  program  rates 
selectable  were  75,  90,  150,  280,  and  450#  per  hour  but  the  actual  program 
rates  were  calculated  for  each  run  from  the  temperature-time  data. 

3.  TEMPERATURE  MEASUREMENT 

The  measurement  of  the  actual  temperature  of  a  material  undergoing  weight 
loss  presents  many  difficulties.  In  theory,  an  ideal  method  for  measuring  the 
temperature  would  be  to  surround  a  thermocouple  bead  completely  with  the 
sample  and  measure  the  thermocouple  millivolt  output.  In  practice,  difficulties 
arise  especially  if  complete  loss  of  material  takes  place  during  degradation. 
In  this  case,  the  thermocouple  bead  becomes  n.  '  and  more  exposed  to  the  heat 
source  as  weight  loss  occurs.  Further,  it  may  iot  be  safe  to  assume  a  correct 
temperature  will  be  dete^  ined  even  when  the  sensor  is  completely  surrounded. 
This  will  be  dependent  upon  the  spectral  characteristics  of  the  sample  and  the 
heating  method  employed  (UV  or  IR). 
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Another  difficulty  is  involved  in  the  simultaneous  determination  of  weight 
and  temperature.  Torsion  of  wires  from  the  thermocoiqjle  to  a  stationary  support 
will  alter  the  mass  reading  or  may  cause  noise  in  the  weight  record.  Methods 
differing  in  complexity  have  been  devised  to  overcome  these  defects  but  none  is 
entirely  satisfactory  (References  19  and  28). 

For  this  work,  it  was  decided  to  measure  temperatures  by  placing  a  thermo- 
couple  in  a  thermowell  as  close  to  the  sample  as  possible.  For  several  of  the 
polymers,  a  series  of  calibration  degradations  was  run.  Sample  temperature, 
measured  by  a  thermocouple  in  direct  contact  with  the  polymer,  was  recorded 
and  corrections  to  the  thermowell  temperatures  were  obtained. 

hi  most  cases,  the  temperature  correction,  A  T,  was  of  the  form 

AT  =  C  +  A/S 


A  and  C  are  constant 
/3  is  the  heating  rate 

The  temperature  corrections  were  usually  in  the  range  of  5  to  15°C.  Similar 
lags  have  previously  been  reported  (Reference  24). 


In  the  range  of  temperatures  over  which  weight  loss  occurs,  severe  temper¬ 
ature  lags  may  occur  (Reference  25)  since  large  quantities  of  heat  are  called 
for  during  an  endothermic  process. 

It  is  realized  that  temperatures  measured  in  this  work  are  only  approximate 
but  it  is  felt  that  by  standardizing  conditions  (sample  size,  crucible  and  furnace 
geometry,  etc.)  and  making  the  corrections  described,  errors  are  minimized. 
It  is  hoped  eventually  to  be  able  to  recalculate  this  data  making  corrections  for 
the  lag  during  the  endothermic  weight  loss.  This  must  await  the  relevant  ex¬ 
perimental  data. 
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4.  PROCEDURE 

The  sample  (usually  100  mg)  was  weighed  into  a  small  quartz  crucible  which 
was  then  suspended  in  a  quartz  degradation  tube  by  a  fine  nichrome  chain 
connected  to  the  balance  beam.  A  counterbalance  was  applied  to  the  opposite 
side  of  the  beam  making  sure  that  weights  at  least  equal  to  the  expected  weight 
change  were  suspended  on  the  beam.  After  the  apparatus  had  been  pumped  down 
to  a  pressure  below  0.1  micron  of  mercury,  the  furnace  which  surrounds  the 
degradation  tube  and  the  programmer  were  switched  on.  After  the  weight  change 
had  occurred,  programming  was  continued  until  a  good  final  weight  base  line 
was  recorded.  From  the  thermogram  which  recorded  both  weight  and  temperature 
as  a  function  of  time,  the  rates  of  weight  loss  as  a  function  of  the  instantaneous 
percent  weight  loss  were  computed  (Section  IV).  In  any  case::  in  which  a  steady 
final  weight  line  could  not  be  obtained  (e.g.,  for  some  aromatic  polyamides 
and  polyesters),  rates  were  based  on  the  initial  sample  weigh-  instead  of  on  the 
total  weight  loss. 

Occasionally  slight  initial  weight  losses  were  noticed  due  to  removal  of 
solvents  or  water  from  the  polymers,  hi  these  uases,  data  ere  taken  from  the 
thermogram  after  the  weight  line  was  again  1:  vel. 

In  an  attempt  to  reduce  procedural  errors,  the  experimental  procedure  and 
apparatus  were  standardized  as  far  as  possible  (e.g.,  furnace,  method  of  tem¬ 
perature  measurement,  sample  size,  crucible  geometry). 

Extremely  useful  Information  on  experimental  methods  and  apparatus  used 
in  thermogravimetry  is  given  in  Reference  25. 
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SECTION  IV 

COMPUTER  HANDLING  OF  THERMOGRAVIMETRIC 
ANALYSIS  DATA 

This  section  is  devoted  to  the  treatment  of  TGA  data  using  Friedman’s 
method  which  is  discussed  at  length  in  Section  II.  There  is  some  similarity  to 
the  treatment  of  isothermal  thermogravimetry  data  (Reference  35).  Despite 
the  possible  repetition,  this  section  describes  all  aspects  of  the  handling  of 
programmed  thermogravimetry  data.  The  first  of  this  section  is  concerned  with 
the  determination  of  the  rates  of  weight  loss  and  the  second  portion  describes 
the  evaluation  of  the  parameters  involved  in  Equation  9. 

1.  COMPUTATION  OF  RATE  OF  WEIGHT  LOSS 

hi  Section  III  the  experimental  procedure  is  described.  Figure  2  shows  a 
typical  recorder  trace  from  the  thermobalance.  The  two  curves  represent 
temperature  (measured  by  a  Chromel/Alumel  thermocouple  located  in  a 
thermowell  close  to  the  sample  container)  and  the  sample  weight  (measured 
electronically  by  determining  changes  in  the  resonant  frequency  of  a  transducer 
caused  by  deflection  of  the  balance  beam).  The  pen  excursions  are  linearly 
dependent  upon  temperature  and  weight,  full  scale  deflections  corresponding 
to  500  or  1000°C  and  10  mg  or  100  mg  weight  change. 

Obviously  the  two  pens  cannot  travel  on  the  same  line  perpendicular  to  the 
time  axis.  A  small  correction  has  to  be  made  to  data  read  from  the  same  line 
to  ensure  that  pairs  of  weights  and  temperature  data  represent  conditions  at 
the  same  time. 

Methods  are  available  for  automatically  converting  signals  from  measuring 
equipment  to  digital  form  for  computer  processing,  but  such  methods  were  not 
on  hand  for  this  work.  Thus  it  was  necessary  to  obtain  the  recorder  traces  and 
to  take  data  from  the  two  curves  either  using  a  mechanical  graph  reader  or 
manually,  hi  the  latter  method,  the  chart  was  taped  to  a  board  and  scales 
graduated  in  suitable  increments  (20th  or  32nd  of  an  inch)  taped  to  both  sides  of 
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the  chart  parallel  to  the  time  axis.  A  ruler  was  then  moved  along  the  time  axis 
and  corresponding  weight  and  temperature  data  recorded  at  200  to  300  positions. 
In  order  to  simplify  data  gathering,  elapsed  inches  of  chart  were  recorded  and 
time  regenerated  during  the  machine  calculations  using  the  chart  speed  (in  inches 
per  minute).  The  zero  time  data  is  taken  from  the  curve  at  some  suitable 
position  prior  to  significant  weight  loss. 

The  first  card  identifies  the  deck  of  cards  that  follows.  The  information 
and  the  columns  in  which  it  is  punched  is  shown  in  Table  IE. 

TABLE  ID 

IDENTIFICATION  CARD 


COLUMNS 

CONTENTS 

1-8 

Any  identification,  name,  etc. 

9-12 

Run  number  (e.g.,  04P) 

13  -  14 

Blanks 

15  -  22 

Date 

23  -  24 

Blanks 

25  -  45 

Alphabetic  information,  polymer  name,  etc 

46  -  52 

Chart  speed  in  inches  per  minute 

53  -  54 

Blanks 

55  -  59 

Pen  offset  in  inches 

60  -  62 

Inches  between  first  and  last  data 

63  -  70 

Blanks 

71  -  73 

Temperature  at  initial  time  reading 

74  -  76 

Blanks 

77  -  79 

Temperature  at  final  time  reading 
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Following  the  identification  card  is  the  data  deck.  The  data  cards  are 
punched  as  shown  in  Table  IV.  A  typical  input  card  deck  is  reproduced  in 
Appendix  L 


TABLE  IV 

CONTENTS  OF  DATA  CARDS 


COLUMNS 

CONTENTS 

1-13 

Identification  information,  name,  etc. 

1  14 

An  index,  LBJ,  to  signal  the  end  of  a  deck 

15  -  17 

Blanks 

18  -  22 

Time  in  inches 

23  -  28 

Weight 

29  -  32 

Temperature 

33 

Blank 

34  -  38 

Time  in  inches 

39  -  44 

Weight 

45  -  48 

Temperature 

49 

Blank 

50-54 

Time  in  inches 

55  -  60 

Weight 

61  -  64 

Temperature 

65 

Blank 

66  -  70 

Time  in  inches 

71  -  76 

Weight 

77  -  80 

Temperature 

Using  the  input  data,  the  computer  assembles  a  matrix  of  time,  weight, 
and  temperature  data  which  is  then  scanned  to  find  the  next  weight  after  a  1% 
weight  loss.  A  number  of  pairs  of  weight  and  time  data  either  side  of  this  point 
is  then  fitted  to  a  quadratic  (using  PLSQ least  squares  curve  fit  subroutine),  The 
quadratic  is  then  solved  for  the  time  taken  for  an  exact  1%  weight  loss  to  occur. 
These  calculations  are  repeated  for  all  integral  percent  weight  losses  193  to  99. 
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Time  and  temperature  data  are  fitted  to  a  polynomial  using  the  PLSQ 
subroutine.  A  single  high  order  polynomial  is  used  to  fit  all  the  temperature¬ 
time  data  whereas  a  quadratic  is  used  to  fit  short  sections  of  the  weight-time 
data.  These  two  different  methods  were  chosen  because  a  quadratic  may  easily 
be  solved  whereas  higher  orders  may  not.  Only  substitution  into  the  time- 
temperature  polynomial  is  needed  here  so  high  orders  may  be  used  to  get  a 
better  fit  of  the  data.  The  weight- time  curve  fit,  however,  had  to  be  solved  to 
interpolate  time  values;  as  a  quadratic  fit  was  dictated,  the  best  fit  was  achieved 
by  using  a  relatively  small  number  of  curve  fit  data. 

After  the  input  data  has  been  compared  with  results  calculated  from  the 
fitted  curves,  the  spurious  input  is  replaced  by  fitted  values  and  the  curve  fits 
recalculated.  The  rate  of  weight  loss  is  then  calculated  for  each  percent  weight 
loss  and  the  average  heating  rate  computed.  A  schematic  representation  of  this 
rate  program  is  given  in  Figure  3,  and  the  complete  program  is  reproduced  in 
Appendix  IL 

Finally  the  results  are  printed  out  in  tabular  form  together  with  graphs 
representing  the  variations  of  rate  of  weight  loss  with  percent  weight  loss, 
weight  with  time,  and  temperature  derivative  with  time.  Results  are  also  punched 
on  to  computer  cards  for  further  processing.  A  typical  punched  card  output 
deck  is  shown  in  Appendix  III  and  examples  of  pointout  data  is  given  in 
Appendix  V. 

2.  COMPUTATION  OF  ARRHENIUS  PARAMETERS 

Friedman’s  method  for  calculating  activation  energy,  etc.,  which  is  de¬ 
scribed  fully  in  Section  II,  requires  the  comparison  of  weight-loss  rates  ob¬ 
tained  from  a  number  of  thermogravimetric  analyses  carried  out  at  differing 
rates  of  temperature  increase.  Comparison  Is  made  between  rates  of  weight 
loss  at  identical  extents  of  weight  loss,  and  activation  energy  may  then  be 
calculated  from  the  slope  of  the  curve  of  log  (rate  of  weight  loss)  against 
reciprocal  of  the  absolute  temperature.  The  lower  the  heating  rate,  the  lower 
the  temperature  will  be  for  a  given  weight  loss. 
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Figure  3.  Schematic  Representation  of  Rate  of  Weight  Loss  Calculations 
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The  output  punched  cards  from  the  rate  program  contain  both  rate  of 
weight  loss  and  temperature  for  each  percent  weight  loss  so  the  program  merely 
selects  corresponding  data  from  each  experiment  and  carries  out  a  least  squares 
straight  line  fit  of  the  log  (rate  of  weight  loss)  and  l/T  data.  Both  the  slope  and 
intercept  of  the  best  line  are  computed.  Log  A  F(W)  values  are  then  calculated 
using  an  average  value  for  the  activation  energy.  The  range  of  weight  loss  over 
which  this  average  is  computed  is  usually  chosen  to  omit  very  low  and  high 
conversions.  A  discussion  of  the  effects  of  changes  in  activation  energy  with 
conversion  is  in  Section  IL 

Figure  4  shows  a  schematic  representation  of  Arrhenius  parameter  cal¬ 
culations,  and  in  Appendix  IV  the  complete  Arrhenius  program  is  reproduced. 
Typical  results  from  all  these  calculations  are  presented  and  discussed  in 
Section  V. 
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Figure  4,  Schematic  Representation  of  Arrhenius  Parameter  Calculations 
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SECTION  V 

WEIGHT  LOSS  OF  POLY(TETRAFLUOROETHYLENE) 

To  validate  the  procedure  for  the  determination  of  activation  energy  and 
order  of  reaction  which  is  described  in  Section  II,  results  obtained  by  this 
method  were  compared  with  those  previously  published  for  the  degradation  of 
polyftetrafluoroethylene),  Teflon.  This  particular  compound  was  chosen  for  its 
relative  lack  of  complications  and  for  the  availability  of  published  information. 

Madorsky  and  co-workers  (Reference  29)  measured  the  weight  of  samples 
of  Teflon  maintained  at  fixed  temperatures.  Piets  of  the  rate  of  weight  loss 
against  the  percent  volatilization  were  linear  between  about  20  and  80%  vola¬ 
tilization  showing  the  degradation  to  be  a  first  order  process.  Their  kinetic 
data  could  be  summarized  by  the  following  equation 


k 


i 


4.7  x 


lrtl8  -80,500/RT  -I 
10  e  sec 


k^  =  first  order  rate  constant 
R  =  gas  constant 
T  =  absolute  temperature 

A  mechanism  for  degradation  involving  thermal,  weak  link,  or  end  initiation 
followed  by  unzipping  of  the  free  radicals  produced  was  postulated. 

Wall  and  Michaelsen  (Reference  30)  confirmed  these  observations  but  sug¬ 
gested  that  a  zero  order  dependence  of  the  rate  of  weight  loss  on  sample  weight 
was  observable  below  about  480°C.  They  presented  data  which  showed  that  at 
460*C  under  nitrogen  the  weight  loss  of  Teflon  is  a  linear  function  of  time  up  to 
about  40%  weight  loss. 

Anderson  (Reference  31)  analyzed  Teflon  thermogravimetry  data  by  the 
method  of  Freeman  and  Carroll  (Reference  11)  and  found  the  degradation,  in 
vacuum,  to  be  first  order  between  450  and  550°C  with  an  activation  energy  of 
75  *4  kcal/mole. 
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Reich  and  co-workers  (Reference  13)  analyzed  Teflon  weight-loss  data 
using  the  method  of  Anderson  and  Freeman  (Reference  12)  previously  discussed. 
They  obtained  activation  energies  varying  between  69  and  74  kcal/mole,  the 
average  being  72  kcal/mole. 

Lee  and  co-workers  (Reference  32)  presented  information  obtained  using 
heating  rates  between  300  and  1200"C  per  hour.  It  will  be  shown  here  that  such 
high  heating  rates  are  likely  to  result  in  large  uncertainties  in  th8  temperature 
measurement  with  consequent  curvature  of  activation  energy  plots.  However 
they  quote  E  between  60  and  69  kcal/mole  and  orders  of  reaction  between  0.7 
and  0.85  depending  on  the  method  of  plotting  employed. 

Carroll  and  Manche  (Reference  33)  re-examined  Madorsky’ a  data  and  deter¬ 
mined  the  aotivation  energy  as  a  function  of  the  conversion  and  showed  that 
between  10  and  80%  weight  loss  the  activation  energy  decreases  from  80  to  about 
46  kcal/mole.  The  decrease  of  E  with  increasing  conversion  was  apparent  for 
both  the  programmed  temperature  increase  and  for  the  isothermal  weight  loss 
of  Teflon  in  vacuum.  The  reaction  was  said  to  be  zero  order. 

Section  Q  gives  a  detailed  discussion  of  the  various  techniques  which  have 
been  used  for  calculating  kinetic  parameters  from  weight-loss  data.  Our  con¬ 
clusions  are  that  each  method  has  serious  shortcomings.  We,  therefore, 
consider  that  the  data  of  Madorsky  (Reference  29)  which  was  obtained  iso- 
thermally  should  be  the  most  reliable  for  comparison  with  the  results  of  our 
investigation. 

As  a  check  on  the  present  differential  method  for  the  determination  of  E& 
and  order  of  reaction,  a  series  of  programmed  temperature  increase,  vacuum 
weight  loss,  experiments  on  100  mg  samples  of  Teflon  was  carried  out  using 
heating  rates  between  45  and  450*  por  hour.  The  polymer  used  was  Du  Pont 
Teflon  molding  powder  Composition  6  in  the  form  of  fine  granules. 

The  first  runs  which  were  carried  out  at  the  higher  heating  rates  gave  very 
high  rates  of  weight  loss  (~10%  per  minute).  With  such  high  rates  of  reaction, 
considerable  temperature  lags  might  be  expected  so  measurements  were  also 
made  using  very  low  heating  rates  (below  150*  per  hour). 
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Figure  5  is  a  plot  of  log  (rate  of  weight  loss)  against  the  reciprocal  of  the 
absolute  temperature  for  50%  conversion.  Plots  for  other  conversions  were 
similar  to  this  one.  It  is  evident  that  a  straight  line  cannot  be  drawn  to  represent 
the  data  over  the  whole  temperature  range.  However,  the  runs  carried  out  at 
150®  per  hour  and  lower  heating  rates  do  show  a  linear  dependence  of  log  (rate 
of  weight  loss)  on  l/T.  The  slope  of  the  line  drawn  through  only  these  four 
points  gives  an  activation  energy  of  69.3  kcal/mole.  A  case  could  be  made  for 
considering  only  the  three  lowest  heating  rates.  The  derived  activation  energy 
would  then  be  increased. 

It  is  probable  that  the  curvature  of  the  Arrhenius  plot  when  the  higher 
heating  rates  are  used  is  due  to  the  large  thermal  lags  when  the  rates  of  weight 
loss  are  large.  In  the  hope  of  bringing  these  results  into  line,  an  attempt  will  be 
made  to  correct  for  these  lags  by  direct  sample  measurement. 

The  computer  printout  for  the  rates  of  weight  loss  for  one  of  the  Teflon 
experiments  is  given  in  Appendix  V  and  Appendix  VI  shows  the  rates  of  weight 
loss  at  each  1%  conversion  for  the  four  lowest  heating  rate  runs.  Appendix  VII 
is  the  computer  printout  for  the  activation  energy  calculations  based  upon  all 
the  rates  quoted  in  Appendix  VI.  Figure  o  contains  all  the  50%  data  from 
Appendix  VI  as  well  as  data  obtained  using  higher  heating  rates. 

1.  VARIATION  OF  ACTIVATION  ENERGY  WITH  CONVERSION 

Figure  6  shows  a  plot  of  activation  energy  as  a  function  of  weight  loss,  the 
data  being  taken  from  Appendix  Vli  Between  4  and  99%  weight  loss,  activation 
energy  varies  between  about  62  and  83  kcal/mole.  The  average  value  of 
69.34  koal/mole  between  10  and  80%  weight  loss  has  been  used  to  calculate 
log  A  F(W)  values  used  in  one  of  the  curves  in  Figure  7.  Careful  Inspection 
of  the  aotlvation  energy  data  reveals  an  approximately  constant  value, 
average  =  63.98  kcal/mole,  between  10  and  50%  weight  loss.  At  greater  con¬ 
versions,  Eft  increases  slowly  to  a  maximum  which  is  maintained  between  65 
and  80%  weight  loss. 

K  has  been  reported  that  Teflon  undergoes  &  change  in  physical  properties 
above  50%  weight  loss.  In  Reference  30  it  is  claimed  that  the  polymer  melts 
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Figure  5.  Arrhenius  Plot  for  Degradation  of  Teflon  at  50%  Weight  Loss 
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Figure  7,  Log  A  F(W)  Curves  for  Teflon  Degradation 
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at  this  conversion.  There  may  be  some  correlation  between  this  phenomenon 
and  the  change  in  the  activation  energy  curve  at  the  same  conversion. 

2.  ORDER  OF  REACTION 

Figure  7  also  shows  the  log  A  F(W)  obtained  using  the  activation  energy 
average  for  the  10  to  50%  weight  loss  (i.e.,  63.98  kcal/mole).  The  straight  lines 
drawn  which  represent  the  20  to  80%  weight  loss  data  fairly  well  have  slopes  of 
0.84  and  0.89,  and  only  slight  maxima  in  the  log  A  F(W)  curves  are  apparent  at 
very  low  conversions.  Thus  the  weight-loss  process  obeys  approximately  first 
order  kinetics.  Other  lines  having  differing  slopes  may  be  drawn  through  points 
representing  more  limited  weight  loss  ranges  but  it  would  be  unwise  to  assign 
real  significance  to  such  slopes. 

3.  PRE- EXPONENTIAL  FACTOR 

The  extrapolated  value  of  log  A  F(W)  at  zero  percent  weight  loss  is  a 
measure  of  the  pre-exponential  factor  (A).  Using  the  low  E  value,  log  A  =  15.7, 

a 

with  the  high  E  ,  log  A  =  17.1.  (Note:  Since  weight  losses  used  here  are  in 

d 

percent,  log  100  has  been  subtracted  from  each  intercept.) 

4.  DISCUSSION 

The  data  presented  here  is  in  fair  agreement  with  some  of  the  published 
information.  The  10  to  80%  average  activation  energy  is  about  10  kcal/mole 
lower  than  Madorsky’s  value.  The  change  in  E  with  conversion  does  not  agree 

a 

with  the  drop  from  80  to  46  kcal/mole  calculated  by  Carroll  and  Manche 
(Reference  33).  It  is,  however,  more  reasonable  to  expect  the  increase  in  E 

Cl 

with  conversion,  as  we  find,  then  to  expect  a  large  drop  in  E  . 

a 

If  the  degradation  involves  unzipping  with  a  long  kinetic  chain  length 
throughout  the  total  weight  loss,  a  first  order  rate  dependency  would  be  indi¬ 
cated,  and  no  changes  in  molecular  weight  of  the  residue  would  occur.  It  has  been 
shown  (Reference  34)  that,  at  500*C,  tetrafluoroethylene  is  the  major  volatile 
product  of  degradation  (95$  but  small  amounts  of  CF^  and  C^Fg  are  also 
produced.  At  higher  temper  .,  the  yield  of  tetrafluoroethylene  is  reduced, 
other  products  being  produced  by  '  "mination  of  short  kinetic  chain  length 
unzipping  processes.  If  the  kinetic  chain  length  is  shorter  than  the  degree 
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of  polymerization,  a  change  in  molecular  weight  of  the  residue  would  take  place 
during  weight  loss  with  consequent  complication  of  the  degradation  mechanism. 

B  is  felt  that  the  method  described  here  for  the  determination  of  kinetic 
parameters  involved  in  thermogravimetry  gives  adequate  agreement  with 
literature  data  for  Teflon  to  Justify  its  application  to  other  systems.  The 
kinetics  of  degradation  of  polymers  which  obey  more  complicated  laws  are 
discussed  in  Section  VL 


SECTION  VI 

DEGRADATION  OF  OTHER  POLYMERS 

In  this  section,  representative  examples  of  results  obtained  during  thermal 
degradation  of  polycondensates  are  discussed.  These  examples  have  been  chosen 
to  show  several  types  of  log  A  F(W)  curves  which  are  derived  from  weight-loss 
data  using  the  computational  methods  described  in  Section  IV. 

1.  POLY  (1,4-  PHENYLENE  SEBACATE) 

Figure  8  shows  how  the  rate  of  weight  loss  for  this  polymer  varies  with  the 
extent  of  conversion.  That  the  mechanism  of  the  degradation  of  this  polymer  is 
more  oomplex  than  that  for  poly  (tetrafluoroe thy lene)  is  shown  by  the  fact  that 
the  curve  of  rate  of  weight  loss  against  percent  weight  loss  exhibits  two  distinct 
maxima,  one  at  about  45%  and  the  second  at  about  90%  of  the  overall  weight  loss. 
Sinoe  separation  of  the  maxima  is  apparent,  the  activation  energies  of  the  in¬ 
dividual  processes  must  differ  appreciably.  The  greater  the  energy  difference, 
the  better  will  be  the  resolution  of  the  rates  of  each  process. 
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Figure  9  shows  that  the  activation  energy  does  indeed  vary  during  the 
degradation  of  poly(l,4-phenylene  sebacate).  Jii  the  early  stages  of  degradation 
(Line  A),  E„  is  about  30  kcal/mole  and  a  gradual  rise  takes  place  until  En  is 

a  a 

dose  to  60  kcal/mole  (hiring  the  final  20%  weight  loss  (Line  B). 

The  curve  of  log  A  F(W)  against  log  (residue  weight-percent),  produced 
when  the  overall  average  E  of  36.1  kcal/mole  is  used,  is  shown  as  Curve  in 

A 

in  Figure  10.  A  change  in  the  slope  of  this  curve  is  apparent  at  conversions 
exceeding  about  70%  showing  the  change  in  mechanism  brought  about  by  the 
commencement  of  the  second  reaction.  As  explained  in  Section  II,  an  erroneously 
high  slope  would  be  derived  from  this  curve  since  this  type  of  plot  should  be 
based  on  the  weight  loss  during  a  single  component  reaction,  here  the  weight-loss 
data  is  based  on  the  sum  of  the  two  component  reactions. 

In  order  to  separate  the  contributions  due  to  each  of  these  two  processes, 
it  is  necessary  to  go  back  to  the  original  curve  of  rate  of  weight  loss  as  a 
function  of  overall  percent  conversion  (Figure  8)  or  the  curve  of  rate  of  weight 
loss  against  temperature  (Figure  11).  By  careful  inspection  and  judicious  use  of 
curve  drawing  techniques,  it  is  possible  to  resolve  the  two  peaks,  from  either 
curve,  into  the  pure  components.  Overlap  of  the  two  reactions  occurs  between 
about  45  and  80%  of  the  total  weight  loss.  When  a  separation  has  been  made, 
calculations  of  log  A  F(W)  for  each  component  may  be  made  using  the  relevant 
conversions.  The  required  activation  energies  are  found  from  the  approximately 
linear  parts  of  the  E  against  weight-loss  curve  (Lines  A  and  B  in  Figure  9). 

Figure  10  shows  the  results  of  such  a  resolution  of  a  complex  weight-loss 
process.  Curve  I  is  the  log  A  F(W)  plot  for  Component  Reaction  L  This  process 
is  probably  random,  the  low  conversion  rise  in  this  curve  probably  being  due 
to  an  early  low  activation  energy  weight-loss  process.  Curve  II  represents 
Component  Reaction  II.  The  drop  in  the  curve  at  low  conversion  may  not  be 
significant  since  this  is  the  region  of  maximum  overlap  with  Reaction  L  The 
slope  of  the  curve  at  higher  conversions  shows  the  reaction  obeys  either  first 
order  or  random  kinetics. 
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Figure  10.  Log  A  F(WJ  Curvet-  Overall  and  Component  Reactions  in  the 
Degradation  of  Pr  v.  i.4-phenylene  sebacate) 
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Figure  11.  Separation  of  Rate  of  Weight  Loss/Temperature  Curves 
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Tims  it  can  be  seen  that  it  is  possible  to  separate  to  some  extent  the  compo¬ 
nents  of  a  complex  weight-loss  process.  However,  where  overlap  occurs  there 
is  some  doubt  about  the  rates  of  the  individual  weight  losses.  This  doubt  is 
reflected  in  the  inability  to  assign  real  significance  to  the  part  of  the  log  A  F(W) 
curve  which  includes  considerable  overlap  of  reactions. 

A  case  could  be  made  from  the  data  shown  in  Figure  8  that  another  compo¬ 
nent  may  be  present  in  the  range  60  to  70%  weight  loss  but  even  if  real  it  would 
be  virtually  impossible  to  achieve  its  resolution. 

A  similar  treatment  of  the  same  data  is  given  in  Reference  23. 

2.  POLY (HE XA METHYLENE  SEBACAMIDE),  NYLON  6.10 

Unlike  many  aromatic  polyamides  (Reference  40), 1  nylon  6.10  degrades 
completely  leaving  no  appreciable  residue.  The  maximum  rate  of  weight  loss 
occurs  at  about  60%  weight  loss. 

Figure  12  shows  the  variation  of  activation  energy  with  the  extent  of  weight 
loss  determined  from  programmed  thermogravimetry  data.  The  acti7ation 
energy  rises  rapidly  during  the  first  10%  weight  loss  and  then  remains  in  the 
region  of  57  kcal/mole  for  the  remainder  of  the  weight-loss  prr  jss.  The  early 
rise  in  E&  can  probably  be  attributed  to  the  removal  of  absorbed  water  from  the 
polymer  or  distillation  of  low  molecular  weight  volatiles.  Using  an  activation 
energy  of  57.2  kcal/mole,  the  log  A  F(W)  curve  shown  in  Figure  13  has  been 
constructed.  The  20  to  90%  weight-loss  data  is  represented  by  a  good  straight 
line  having  a  slope  of  1.02  indicative  of  a  random  or  first  order  decomposition 
mechanism.  The  downward  curvature  of  the  line  at  low  conversions  tends  to 
indicate  a  random  process  is  operative  especially  as  the  maximum  occurs  close 
to  25%  weight  loss  (Reference  5).  Other  workers  have  concluded  that  the  same 
mechanism  describes  the  degradation  of  other  polyamides  but  a  possible  ionic 
hydrolysis  process  may  occur  simultaneously  (Reference  38). 

Several  other  aliphatic  polymers  have  been  studied  and  found  to  yield 
similar  log  A  F(W)  curves  end  mechanistic  interpretations.  The  degradation  ul 
several  aliphatic  polyesters  and  polyamides  is  discussed  in  References  36  and  39, 


Figure  12,  Activation  Energy  for  Weight  Loss  of  Nylon  6. 
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SECTION  vn 
CONCLUSIONS 


The  method  of  Friedman  has  been  shown  to  be  applicable  to  a  wide  variety 
of  polymer  degradations.  It  can  detect  a  single  first  order  weight  loan,  one 
where  random  decomposition  is  suspected  (as  noted  by  a  maximum  in  the  rate 
law  curve)  and  can  provide  some  resolution  of  a  weight-loss  curve  showing 
more  than  one  rate  maximum,  ft  is  our  conclusion  that  only  through  the  use  of 
several  thermograms  at  different  heating  rates  can  this  information  be  obtained 
reliably. 
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C  PROGRAM  TO  DETERMINE  RATES  OP  WEIGHT  I155T  ATEjnE  PERCENT  vcElGHT  LOSS 
C  INTERVALS  FROM  THERMOGR AV I METP I C  ANALYSIS  DATA. 

C  INPUT  TEMPERATURES  ARE  FITTED  TO  A  FIFTH  DEGREE  POLYNOMIAL  USING  A 

C  LEAST  SCUARES  SUBROUTINE  (PLSO). 

c  heights  corresponding  to  short  temperature  ranges  are  fitted  to  a 

c  cuauratic  by  plsq. 

T - TfsTPUT^ wTT TCFTS  “DITFERTING- PRDM  FITTED  XT N E~ RT M 0 R E"T H A R  UNIT  P  CPC EnT0T  Th E 

C  TOTAL  HEIGHT  LOSS  ARE  REPLACED  BY  THE  CURVE  FIT  VALUE. 

C  OUTPUT  DATA  IS  PUNCHED  ON  TO  CAROS  FOR  FURTHER  PROCESSING  (TO  CALCULATE 

C  ACTIVATION  ENERGY  ETC). 

C  INPUT  DATA.  SPCED=ChART  SPDED  TFT  INCHES  PER  MINUTE,  OFF SET* PEN  SEPARATION 
C  IN  INCHES.  ITI-INITIAL  TEMPERATURE  READING,  I TF=F  INAL  TEMP.  READING. 

C  TNDEX  LBJ  IS  SfT  EW5r'TD-T~TTrTI0¥AL  LAST  "C*KD  DF  ATTOiT  . 

C 

c 

DIMENSION  T 1 1500) ,H( 1500) ,TEMP( 1500  > , WWI 1 05 ) , TNW( 105 ) , DWOT ( 105 ) , 

•  X(5C),v(5O),C(l0),A(l l.PLOT' 50, 105) ,Z<  50 V.TDER  t ID5),TP0LY< 105) , 

.RTEMP( 105) 

1  READ  (5,1000)  I O,  DATE  1,0  ATE  2  ,  COM  1 ,  COM2  , CQMJTCOMa*  SPEEDTOFFSET, 
•XRANGE,!TI ,ITF 
CHART  =  2.5A5A5A*SPEED 
WRITE  (6,3000) 

WRITE  (6,3050)  I D, D AT E 1, 0 AT E 2, COM  1, COM 2, COM 3, COMA 

WRITE  (6,3060)  SPEED 

WRITE  (6,30651  CHART 

WRITE  (6,3070)  OFFSET 

WRITE  (6,3080)  XRANGE 

WRITE  (6,3090)  I T I 

WRITt  (6,3100)  ITF 

N  *  0 


1 

5 

6 
7 
B 

9 

10 
11 
12 


START  LOOP  TO  READ  IN  DATA  CAROS 


DO  30  1=1, 2995, A 
N  *  N  + A 

READ  (5.1010)  LBJ,T( I) , W( I !, TEMPI  I) .  T( I ♦ 1 ) , Wl I  ♦  l ) , TEMPI  I ♦ 1 ) , 

'.TTIV2 )  ,wll*2 )  ,TfHPTI*2T .TITdl , w (I ♦  3 ) . Tfc  1*31  " 

t  «  time  (in  inches),  w  «  weight,  temp  «  temperature 

17 

IF  (L8J.E0.1I  GO  TO  AO 

30  CONTINUE  _ 

AO  NL  *  N-A  . 

C 

t  START  LOOP  TO  CALCULATE  NUMBER  OF  DATA  READ  IN 

C 

DO  50  I *NL  , N 

IF  ( T ( I I.NE.O.O)  GO  TO  50 
N  »  1-1 
GO  TC  51 

50  CONTINUE 

51  JJ  -  . 030*F LPAT ( N ) 

LL  ■  MAXO(JJ.IO) 

C 

C  JJ  »T  PERCENT  OF  no.  cf  DATA  SETS  rEaO  IN 


50 


ooo  ooon  ooo  ooo  ooo  ooooo 


\  FAIL- TP.-'*- 1*1 
r’sri  II 


TGA 

PLOT  -  EFN  SOURCE  STATEMENT  -  IFN1S)  - 

C  LL  =  NO.  OF  CURVE  FIT  POINTS  (LAtCR  *  KM) 

C 

WRITE  ( c , 3170 )  LL 
WRITE  (6,3C10)  N 
NN  =  LL 

K  =  POLYNOMIAL  ORDER, NEEDED  FOR  PLSU  SUBROUTINE.  LIST  * 
ANALYSIS  OF  PLSC 
D  *  TOTAL  WEIGHT  LOSS 

K  =  2 
LIST  *  0 

LINDA  =  l  __  _  _  _ _ _ 

b  =  win  -  w(n> 

00  i 5  1=1, N 

W(I)  =  100. -1100. *(W( 1 l-KINII/D) 

CONVERT  TIME  IN  INCHES  TO  MINUTES 

Till  =  T ( I  I/SPEED 
55  CONTINUE 

DT  =  PEN  SEPARATION  IN  MINUTES. 

DT  =  SPEED*OFFSET 

START  MAJOR  LOOP 

DO  ICO  NW  *  1,99 
58  II  =  LINDA-1 

WW(NW)  =  FLOAT(NW)  _ 

SCAN  WEIGHT  DATA  FOR  ONE  CLOSE  TO  BUT  JUST  GREATER  THAN 
LOSS.  II  = I  NOE  X  OF  THAT  POINT 

DC  tO  I =L I NDA , N 

II  =  11*1 

IF  ( W ( I l.GT.ww(NW) )  GO  TO  70 
60  CONTINUE 
70  LINDA  =  I  I- (LL/2) 

LINDA  =,  INCEX  OF  FIRST  DATA  TU  BE  USED  BY  PISQ 


DO  eC  J= 1 , L L 
JI  =  L I  NO  A ♦ J- 1 
X  (  J  »  =  T( JI  ) 

Y  (  J  1  *  W  (  J  I  ) 

7! J)  *  TEMP( JI  ) 

8C  CONTINUE 

c 

C  CURVE  FIT  CF  TIME  AND  WEIGHT  DATA 

C 

CALL  PvSW  (X,Y,NN,K,C,LIST,EMAX, ERMS,EMEQ) 
K  K  =  l 

c  __  _ 

C  START  LOOP  TO  ChECK  FUR  RAD  INPUT  DATA 


10/30/63 


0  FOR  NO  ERROR 


ONE  °ERCENT  WEIGHT 
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EFN  SOURCE  STATEMENT  -  IFN(S> 


10/30/68 


C 

- jnrin~j*TTLT 

JI  «  UNOA+J-1 _ _ 

C  WE  »  WEIGHT  CALCULATED  FROM  POLYNOMIAL 

c  '  . . . . 

WE  *  Cm*TIJll«2*CI2l*TI JI)+C(3) 

~ C - - - -  - - 

C  COMPARE  CALCULATED  AND  ORIGINAL  DATA 

IF  (ABSIWE-WUm.GT.l.:  GO  TO  82 

- - - CO  TC  101 - -  -  -  ' 

82  WRITE  (6,4000)  JI,M(JI).Wg _ _ 

C  REPLACE  BAC  OATA  BY  CALCULATED  VALUES 

~  C  ~  '■  '  "  "  '  ~ .  '  . .  "  . . . .  '  . 114 

_ W(  JI)  *  HE _ _ 

81  CONTINUE 

- 50" TC  "(83,58fVKK  - - 

c 

~  '  CHECK  FOR  IMAGINARY  ROCTS  IN  SOLUTION  OF  QUADRATIC 

C 

83  SCREW  *C  <2 ) *C ( 2  )“4 .  0*C  (  l )  *(C I  3 ) -WW  ( NW ) ) 

IF  (SCREW. LT. 0.0)  GO  TC  90 _  _ 

C  USE  REAL  RCOT  TO  DETERMINE  TIME  CORRESPONDING  TO  EACH  PERCENT  WEIGHT  LOSS 
C 

TNWINW)  *  (SQRT(C(2)*C12)-4.0*C(1>*(CI3)-WWINW) ))-C(2> )/(2.0*C(l> ) 

c  . 

C  DWDT  *  RATE  OF  WEIGHT  LOSS 

~X  "  '  -----  - . '  . .  126 

DWDT ( NW )  «  2.0*C(1»*TNH(NW)  ♦  C(2) 

GO  TC  100  "  '  . . 

90  TNWINW)  *  T(II-l) 

OWOT(NW)  •  0.0 

TA  »  Till) ASPEED _ 

C  WRITE  OUT  IDENTIFICATION  AND  LOCATION  OF  BAD  DATA 

C  ”  . . . 

WRITE  (6,3160)  NW, 1 1 ,TA,W( 1 1 )  136 

'  ICO  CONTINUE 
K  «  5 

~Z - - *  '  ■' 

C  CURVE  FIT  CF  TIME  AND  TEMPERATURE  DATA 

'  c  '  " . . . ' . . . . . 


CALL  PLSQU,  TEMP,  N,K,C,  LIST,  EMAX, ERMS, EMEQ >  142 
WRITE  (6,3100)  EMAX  140 
WRITE  (6,5200)  ERMS  144 

~wrtte  rerrsoD  ~  ' . "  ■  _ . "  ~  145 

WRITE  (6,5400)  (C(I), 1-1,6)  146 
WRITE  (6,3110)  153 
00  120  NW* 1,99 

c  . . 


C  CT  *  CORRECTED  TIME,  USED  TO  FIND  TIME  AT  EACH  TEMPERATURE 
C 


52 
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10/30/68 

CT  *  TNW( Nh )-DT 

c 

TSTOR  *  Ctl) 

c 

c 

LOOP  TO  EVALUATE  TEMPERATURE  POLYNOMIAL  FOR  EACH  VALUE  OF  CT 

200 

00  200  1*2,6 

TSTOR  *  TSTOR*CT+C ( I ) 

TPOLY(NW)  *  TSTOR 

TSTOft  *  5. *C ( l ) 

00  250  1*2,5 

J  *  6-1 

C 

250 

TSTOR  *  TSTOR*CT+FLOAT(J)*C( I ) 

~T 

C 

TDER  *  TEMPERATURE  DERIVATIVE 

RTEMP  *  RECIPROCAL  OF  ABSOLUTE  TEMPERATURE 

C 

TDER(NW)  *  TSTOR 

RTEMP(NW)  *  1. 0/ ( T POL Y(NW 1+273. 16) 

WRITE  (6,3120)  NW.OWOT (NW ), TPOLY( NM ), TDER ( NW ), RTEMP (NW) 

172 

120 

CONTINUE 

STDER  *  0.0 

C 

C 

CALCULATE  AVERAGE  TEMPERATURE  DERIVATIVE  (AVE) 

C 

DO  125  1*1,99 

125 

STDER  «  STOER  ♦  TDER(I) 

CONTINUE 

AVE  a  STDER/99.0 

WRITE  (6,3125)  AVE 

C 

C 

SET  UP  DUMMY  POINTS  FOR  GRAPH  PLOTTING  SUBROUTINE  (GP) 

C 

WW(IOO)  *  0.0 

187 

DWDT(IOO)  *  0.0 

TNW(IOO)  *  TNW ( 99 ) 

TOER(IOO)  *  0.0 
!  WW(lOl)  -  100.0 

|  DWDT(lOl)  *  0.0 

!  TNW(lOl)  -  TNW ( 99 ) 

1 

TDER ( 101 )  -  TDER(99) 

WRITE  (6,3000) 

188 

((RITE  (6,3130)  ID 

L  «  3 

189 

|  IS  -  5 

}  . _ LW  »  101 

!  LN  «  50 

! _ H  »  101 _ 

OATA  A/1H. / 
JN  -  1 _ 


PLOT  GRAPH  OF  RATE  OF  WEIGHT  LOSS  AGAINST  PERCENT  WEIGHT  LOSS 

CALL  GP  (WW.OWOT ,L,LS,H, JN. LW, LN, A. PLOT  I _ 

WRITE  (6,3000) 

WRITE  (6,3140)  ID _ 


C 


196 

197 


AFML-TR-F8-181 
Part  0 


Tga  '  ’  10730/68 

PLOT  -  EFN  SOURCE  STATEMENT  __•*  I FN ( SI  - 

C  PLOT  GRAPH  OF  PERCENT  WEIGHT  LOSS  AGAINST  TIME 

X -  - - -  lw 

CALL  GP  (TNW.WW,L,LS,M,JN,LW,LN,A,PLOT)  199 

WRITE  (6,3000)  100 

WRITE  (6,3150)  10 

C 

C  PLOT  GRAPH  OF  TEMPERATURE  DERIVATIVE  AGAINST  TIME 

X - 201 

CALL  GP  (TNW,TOER,L,LS,M, JN,LW,LN,A,PLOT)  _  202 

OWDT(IOO)  *  0.3 
TPOLY ( 100)  -  0.0 

C  "  ~  . "  .  " . ■ . "  ~  "  "  “  " 


C  PUNCH  OUTPUT  CAROS  CONTAINING  PERCENT  WT-  LGSS(NW)  THEN  THREE  PAIRS  OF 

"X - TEKPERTIWr  SND~FT5T T  OP  W XT GHTTO $X~XSTA 

C 

00  150  NW= 1 , 100,3 

PUNCH  5000 ,10, NW ,DW0T ( NW ) , TPOLY ( NW ) , DWDT ( NW*1 ) , TPOLY (NW+1 ) , 
.0WOT(NW+2T,TPOLY(Nw+2)  206 

150  CONTINUE 

t~ 

C  LOOK  FOR  FURTHER  SETS  CF  DATA  IF  NEXT  CARO  CONTAINS  A  ONE  IN  COLUMN  10 


READ  (5,1020!  MORE  215 

IF  (WORE. EC. 1)  GO  TO  1 

_ STOP _ _ _ _ 

1000  FORMAT (8X,A4,2X,A6,A2,2X,A6,A6,A6«A3«F7.4,2X,F$.4,2X,F6.2,3X, 

.I3.3X.I3) 

1010  F0RMAT(13X,I1,2X,3(F6.2,F5.1, 1X,F4.0)»F6.2,F5. 1.1X.F4.0) 

1020  FORMAT  (9X.I1) 

1030  FORMAT  (121 

.3000  FORMAT  (1H1) _  _  _ 

3010  FORMAT  ( 10X , 21Hf OTAL  NC  OF  POINTS'"-  ,!4f 
3050  FORMAT (5X,A4,BX,A6,A2,13X,A6,A6,A6,A3//) 

3060  FORMAT  (10X.33HCHART  SPEED  (INCHES  PER MINUTE)  -,F6.4) 

3065  FORMAT  ( 10X , 38HCHART  SPEED  (FURLONGS  PER  FORTNIGHT)  «,F7.4> 

3070  FORMAT  (IOX.IThOFFSET  (INCHES)  * , F7 .«!  . 

3080  FORMAT  ( 13X, 17HXRANGE  (INCHES)  »,F7.2) 

3OWT0RM T  TIOx. I BhTNTTI At  TEMP  ( C )  -,I5T 
3100  FORMAT  ( 10X  ,  16HF I NAL  TEMP  (C)  *,I5) 

3106  FORMAT  ( Tf X , 1 5 « 4x , F  9 .3, 4  X , FT. 2 » 4x7F 6721 

3110  FORMAT  ( //3X, 11HWE IGHT  LOSS , 6X, 8HDWDT ( NW ) , 13X , 5HTP0LY ,6X, 

.4HTDER,llX,5HRTtMP) 

3120  FORMAT  <6X ,  1 3 , 10X , El 2 . 5, 7X , F9.3, 2E 1 5 .5 ) 

3123  FORMAT  (  V10X.27H  AvERAGE_7XMP~15FRrvSTIvr'»_7FT5;51 
3130  FORMAT  (10X.19H0WDT  VS  WEIGHT  L0SS,20X,A4) 

3149  FORMAT  "TTOX  « 19HWE IOhT  LOSS' VS  7ImE»26X.A4> 

3150  FORMAT  ( 10X , 12HTDER  VS  T I  ME , 20X , A4 ) 

3160  formaT(2x, IThSCrEW  LESS  Than  0 ,T5x » 3hnw»,  13^ 1 OX ,  3Hl I« ,1 4 , 15  x , 

•2HT-.F6.2, 10X.2HW- , F5 . 1 ) 

HTX'FOR'MAT  M5x'i?3HNTrXFTT3_  IN  CURVE  TTT“-~Tm 

4000  FORMAT  (10X.9HAT  PT  NO  ,I4,10H  WEIGHT  ■  ,F5.1,13H  REPLACED  BY  , 

.F5.1)  ' 

5000  FORMAT  (2X.A4, I4,E13.5,F6.1,E13.5,F6.1,E13.5,F6.1 ) 

5160  FORMAT  (lOX.lTHMAX  fEMP  ERROR  »  .FI6.61 

5200  FORMAT  (10X,30HTEMP  ROOT  MEAN  SQUARE  ERROR  «  ,F10.6) 

T366'T0RMATTIOX,15HfEMP  POLY  COEfF)  _ 

TGA 

PLOT  -  EFN  SOURCE  STATEMENT  -  JFNISI  -  _  _ _ _ 

5400  FORMAT(13X,F10.6)  _ 

ENO 
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4  0.13932E 

7_ 0.2244 6 t 

13  0.37753c 

_16  0.44440c 

19  0. 53156 E 

22  0.57782E 

25  0.61520E 

TB — 3TS1T4T5F 
31  C.74898E 

34  0.76522c 

37  G.£T0c'4Cc 
40  0.85006c 

43  0.8746SE 

46  0  •  6  H  / 4 3 1 

49  0.94210E 

52  C.97623E 

*5 5  * ’0  .’9 7  1 7 8  l 
58  C.1C131E 

61  C. 10400c 
T5 — uTKnrror 

67  0.10344E 

70  U.1C387E 

"T3  "0. to 3  21 E 
76  C.1CCUE 
( 5  w.IOC^jc 
*82  3V907Z4L 

85  0.86745c 

86  0.79915E 

"91  0.72T52E 

94  C.59624E 

97  C.41373E 


TEFLON  DEC 

-0D~ZZ^7T 
01  529.2 
01  536.7 
ii  545.fi 
01  5  46.5 
01  554.0 

HnT[.T 

01  560.5 
01  563.3 
0 1  566.0 
01  568.3 
01  570.4 

~0T  "5  7  2Y6 
01  574.6 
01  276.-; 

01  580.2 
01  561.9 

'01'  '582.7 
02  565.4 
02  58710 
02  I oo.u 
02  590.2 
02  591.8 
02  ~5’3’3V4 
02  595.1 
02  556.8 
01  3  rj  8.5 
0  i  6  00.5 

01 _ o  02 . 4 

01  604.6 
01  607.2 
01  610.7 


,A7 1 CN _ 

.7 6  4~67  c 
>17 160c 


.  j 5 2 6oL 
.38G77E 
.63133c 

mmrr 

.73  789E 
.77316E 
.62437c 


TO  j  1 6  •  1 ' 
01  552.5 
01  ‘541.1 

HTswrr 
01  351.3 
0 1  _5  55 . 2 
0  I  -•  C  .  6 
01  361.3 
01  564.2 
'0I“366. r 
01  569.0 
01  571.2 
'01  373.3' 


Tncm 


j68w  0*  6  i*  7  < 


'.7  '  O'. 

.4  0. 


G  J V J JL 

Oo  762  3c 
0.95691c 
0 . 1020P.E 
0. 10414c 
o .  io2o4c 
0 .103^8- 
0.10423E 
0.10237c 
0.10x516 
0.99414c 
0 . 5  h  0  1  0  c 
0 .547o6c 
0.7737^c 
'Oi  6  c  TBIT 

0 .34449c 
0.3328CE 


0  *  5  c  0  .  u 
01  562.5 
01  584.2 
02  585.9 
02  567.5 
T2  385.1 
02  590. c 
02  332.4 
02  '594.0 
02  395.7 
0 1  5  9  7.3 
CIT  3  3  5 .  CT 
01  601*0 
01  6  0  ■>  •  1 
Cl  5"0 1". 4 
01  603 .2 
01  612.3 


*Oc  01 

ri  oi 
luTTT 
jiE  01 


9l"0  1 
1 9c  01 


5  9  E  01 
>3E  01 
V3E  Cl 
j 6c  Cl 
3  3  c  01 


5  7  4  4  s  : 
10038c 
10281c 
10376c 
1 033  Z"l 
10331: 
10435: 
lCOvOc 
10124: 
3 79371 

~ 1 0  3  37>~ 

6  c  3c0: 
75096c 
"64  573": 
46352c 
23131C 
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Y(J4  “  06/13/6  9 

PLOT  -  EFN  SOURCE  STATEMENT  -  IFN(S)  - _ _ 


" C - PVDfflAMtfE  TF  0  E  T  E  E  M IN  r  TCA  'PArAHE TEff S'ffV  FRIEDMANS  METHCP 

C  PROGRAMME  ACCEPTS  OATA  CAPOS  HAVING  THREE  SETS  OF  DATA  PER  CARD. 

“C - L  AST  C«D"0F  ?,ACH"TC{rK,'MTJ5T  HAVE  A  ONE  "IN"  COLUMN I.  LAST  CARD  OF 

C  LAST  OECK  FOR  ONE  °OLVMER  SYSTEM  MUST  HAVE  A  TWO  IN  COLUMN  1  INSTEAD 

T - TO  'ffWT  SECOND  S F T”  OF"  "DEC  TC5  V " BUNCH  T  CARD  WITH  A  THREE  IN  COLUMTTT 

C  AND  PLACE  BETWEEN  SETS 

T - AT'  THE . FTO'ITF  ALL  ""DeCK  5"  PU'SCF"~A"  BLANK  CARO — THEN  AN — TCTJF 

C 

~C - SVHWLS — RffDT- A'TTiyT'-^T5HT"'TnT3nirr"PrT ?  M  B  =  &  E  C  rWCTATTlT'  AB  SOTUTE 

C  TEMPERATURE,  RATE  *  LOG  RATE  OF  WEIGHT  LOSS.  SLOPE  =  SLOPE  OF  ARRHENIUS 

T - PT0T7'""PRnF X "  ■  PP'E- ETPDN E NTTA'L" "FA  CTOTPT~P LOT  =  DIMENSION  OF  OP  SUBROUTINE 

C  ACTE  *  ACTIVATION  ENERGY,  X  AND  Y  REPRESENT  DATA  TREATED  BY  GP 

"C . TPOrr  w ■  INPUT  ■  TT: M’P'ETrATIJPTSir r^-TDcNTT'FTC3TT'D'N7“A- =~NTf  “OF  "SYMSOLT~TfriTP 

C  A A  *  PERCENT  WEIGHT  LOSS,  AFW  »  FUNCTION  FROM  FRIEDMANS  EQUATION 

~C - Fw~WPR  AGf-AFW,  A  A'  LO'GTPc'ft'CTNT  BE  5'IW  TV"  wF"  =  AWPA'CE'TFW - 

C 

- BWKiSTtiN  ■'DftDTn'FF.'lCT.U  TEMPI  !  30',  191 ,  R  2  T?  f  !"5S '» V3TVSLFPF  ( 1W1 , - 

.PR  EX (100), PLOT (50»100),ACTE(100),X(10)»Y(  10), 

“■  “ “ iTPDLVn DDV1 5T,T DCTFTTA  (!) , AA ( 100) ,  AFwTl  ?5) ,  FWT  1  Fil  ,'B (H C3 ST  i¥Fi '35') , 

.SPS(tOO)  ,SDS<  100), SDK  100) 

1  READ  ( 5, IcOO )  IG,CnMl,C0M2, COM*, COM4, C0MS,Cf)Me  ,ComT,Ccm8 
WRITE  (6,3000) 

WRITE  (6,11001  IG,  COM)  .COM2,  COM3,  COM4 , £3*5 , C0m6 ,  COM?,  CC«B 

2  J  -  0 
"lo  J  «  J+T 

C 

~C  StapT  lOop  To  read  In  OAta 
c 

~  00  20  NW  a  l,  •rf7? 


T 


LBJ  =«  V  TN~CnLUMNT~OF  LAST-FARO  OF'  A  DTCk  ,  '  L  AST  CARD  OF  LAST  irETTK'To'?? 
ONE  POLYMER  SYSTEM  NEEOS  LRJ  =  2. 


R  E AO  I*. 1200  I  LRJ,ID( J ), I W,DWOT(NW, J ) , TPOL Y(NW, J ) ,DWDT (NW+1 , J) , 
.fPOLYl NW+1 ,J),0WDT(NW  +  2, J) ,TpOLY(  NW*2, J) 


CHECK  THAT'  INPUT  "CARD'S  AB  ="TN'CONSlC0T  t  V?  OF  Of  R 

10 

(  F  (IW-MWn,4,3 

3  WRITF  (6,1900)  NW, I0( J ) , IW  21 

STOP . 

4  AA(NW)  =  FLOAT(NW) 

A  A  (  NW+1 )  a  FLO  AT(  MW  +  1  )" .  '  . . . .  ~  ' 

A A( NW4? )  =  FLOAT ( NW+2 ) 

IF  (LBJ.'EQ.l  TOO  TO  ir 

20  IF  (LBJ.FQ.2)_G0 JO  25 _ 

?5X J  *J  . .  ” 


WOITg  LIST  OF  sun  IDS 

w rite  T 4 , f Mo oTTToT II, 1  =  1, j) 

C  HE  CK  S  0  R~  A  T  T  c  AS  t  T  HP  EE  OATA  OcO<S 
Ic( J-’)  30,35,36 


3  6 


Best  Available  C( 


58 


oo  o  ooo  u  u  u  ooo 


AFML-TR-68-181 
Part  H 


TOA 

PLOT  -  EFN 

SOURCE  STATEMENT  -  IFN(S) 

30 

WRITF  (6,2000) 

GO  TO  300 

35 

WRITE  (6,1500) 

TSUM  =  0 

N  =  C 

S  °R  EX  =  0.0 

START  LOOP  TO  CALCULATE:'  LEAST  SQUARES  LINC  nc  LOG  (PATH)  VS.  FT5M? 

DO  45  NW  =  4,08 

SIJMXX  =0  _  _ 

SUMYY  =0 
SUMX  =  0 
SUMY  =  0 
S  UMXY  =  0 
on  40  K  =  1,J 

CHECK  POP  7FprJ  RATCS 

I  F(  DW0T ( NW»K ) .LT . 1 .0E- 1 0)  OQ  TO  *,5 

R  ATF  (NW, K )  =  AL0010(DWPT{MW,K)  )_  __  5-! 

RTEMP(NW,K)  =  1.0/<tpplY(NW,K)*273.16| 

S'JMXX  =  PARTIAL  SUM  np  X  SQUARED  ETC. 

SUMXX  =  SUMXX  ♦  prc^pi NW,K  )**2 
SUMYY  =  SUMYY  ♦  ( c ATE (NW t K ) ) *  * ? 

SUMX  =  SUMX  ♦  OTKMP(SW,K  ) 

SUSY  =  SUMY  ♦_RATP(MW,K) 

40  SUMXY  =  SUMXY  ♦  RTEMP(NW,K)*faT-("nw',KJ 
c.o  rn  5  5 

SET  UP  DUMMY  points  FOR  OP  IF  A  0W0T  VALU:  IS  TEP1' 

55  ACTE ( NW )  =  0. 

PRE*(NW)  =  0. 

PATE  (  NW ,  K  )  =  0. 

R  TE  MP  (NW.KI  -  0.0015 
00  TO  45 

55  SLOPE  (  NW  )  =  <XJ*SUM*Y-SUMX*SUMYl/<  XJ'SUMXX-  S'P'X»*2) 

SPS  (\'W  )  =  (  (SUMYY-  (  S  I  lMY  *  SUM  Y/  X  I  )-  (  (  X  J*  SUMXY-  SUM x * SU^Y ) * *? / 

.(  XJ*XJ*SU“XX-XJ*SUMX*SUMX)  ))/( XJ-P.C)  ) 

ALPHA  =  (  SPS(NW)  /  <  SIJMXX-  (  Si(MX*StlMX/XJ  )  )  )  *4.576 
IF( ALPHA)  5S,Sfl,S7 


57 

sns(NW)  = 
00  TO  5P 

SORT ( ALPHA) 

-IS 

5R 

SDS(NW)  = 

".0 

59 

«ETA  =  (SPS(NW)*SUMXX/(XJ»'SUMXX-SUMX*SUMX)  ) 

I F( RETS ) 

62,63,51 

61 

SDI(NW)  = 

00  TP  63 

SOP  T ( RCT  A  ) 

3, 

62 

SOI (NW I  = 

3.3 

53 

ACTE(NW) 

=  -SLOPE(NW)*4.5T6 

PPFX(NW)  =  (  SUMXXOSUMY- SUMX*SJMXY) /(  X  J*SIJM*  X- SUM  Y  2  | 
I F(NW.LT.IO)  on  TO  45 
IF(NW.OT.PC)  00  tp  45 
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TGA  "  Ci>/1?76» 

plot  -  :-:fn  source  statement  -  ifnis)  - 

T  SUM  =  TSUM-SLOPE(NW) 

•spur  « TfPBix-rPTTjnwT  ~ - - - -  - —  • 

N  =  N-M 

45  CONTINUE  ' 

CALCULATE  AVERAGE  ACTTVA'T  TON  ESTEP  (TV  AN  IV  P®E“"EVPr'N'EVTTffL '"F'A’CTUP”" 

- A VPPE <  •»  'SoPFTT^TOATTTn - - -  - - 

AVEA  *  TSUM/FLOATIN) 

A  VAC  T  E  =  AVEA*1.9R7*2.303 

START  LOOr  TO  CALCULATE  AEw 

CO  TO  NW  *  4", p A 
?  =>  0 

00  00  K  =  1,J 

ft  FW ( K  )  *  RATE!  NW*  K  )  A  Vf  A  *R  T  =  MO(  NW  ,K  ) 

90  I  *  7  ♦  AFWlK) 

FW(NW)  •-  I/XJ 
WN  »  FLOAT  fNW ) 
or,  *  A10G1C< 100.-WN) 

SO  1  0 

on  pi  k  =  i,j 

P 3  SO  *  SO  ♦  (FW(NW)-4FW<K) )**? 

VK  *  J-l 

SOAFW  =  SGRTISO/Vk) 

WRITE  OUT‘RESULTS  oFOfENT  WT.  LOSS,  ACTIVATION  rNFFGV,  PP  c -F  *  PON =  'jt  I  At 
FACTOR,  AVERAGE  FW.  AND  ST  AND  aRU  DEVIATIONS,  ALSO  LCr.  WEIGHT  REMAINING!  GG  ) 

70  V R  I  T F  (6.I400I  NW,  'CTE(NW)  ,SOS(NW)  LPP_F  M  NW )  ,  SC  I  <  NW  )  ,  FW  <  NW )  ,  SDA  F  W  , 

".GG 

WRITF  ( 6i 1 425 )  AVACTE 
WRITE  (6,1435)  AVPREX 
WRITc  (o,1440) 

S  FT  UP  INFORMATION  Ft?  r,o  SUMROUTINt,  SF r  othtf  Ofi(T,RA«S 

L  *  1 
L  S  *>  5 
l W  *  l^n 
LN  *  5  n 
M  «  J 

DATA  A/1H./ 

JN  »  1 

START  LOOP  for  BLOTTING  C.RARMS  AT  10  P E DC  -NT  WEIGHT  l  'O  INT'PVAIS 

00  200  NW  *  10,90,10 
00  ICO  X  »  1,J 
X(«)  *  RTF  MR ( NW , X  ) 

100  V(K>  »  T  A  T  F ( NW , K ) 

WRITE  (6,3000) 

WRITE  (*-,l'00)  NW 

C  PLOT  r,  0  A  o  H  CE  LOG  (WATR  OF  WFI'.HT  L  t)  s  C  )  AGAINST  ‘TpJ.viL 
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TQA  06/13/68 

.  PLOT  -  EFN  SOURCE  STATEMENT  -  IFN(S)  - 

C 

OF  TEMPERATURE 

c 

200 

CALL  GP  (X,Y,L,LS,M,JN,LW,LN,A,PLOT) 

161  | 
163 

M  a  100 

WRITE  (6,3000) 

167 

c 

WRITE  (6,3100) 

•  c 
c 

PLOT  GRAPH  OF  ACTIVATION  ENERGY  AGAINST  PERCENT  WT.  LOSS 

168 

CALL  GP  (AA, ACTE,L ,LS. M,JN,LW,LN, A, PLOT) 

WRITE  (6,3000) 

16  S 
17C 

c 

WRITE  ( 6 , 3 2C 0 ) 

c 

c. 

PLOT  GRAPH  OF  PRE-EXPONENTIAL  FACTOR  AGAINST  PERCENT  WEIGHT. LOSS 

171 

CALL  GP  ( AA,PREX,L,tS,M, JN,LW,LN, A, PLOT) 

WPITF  (6,3000) 

172 

173 

WRITE  (6,3300) 

DO  75  I«l,87 

174 

78 

891  I)  «  A10G10( 10C.-AA< 1  +  3 1 ) 

WF(  I)  a  FW( 1  +  3 ) 

175 

LW  *  05 

M  a  87 

C 

c 

PLOT  GRAPH  OF  LOG(AFW)  AGAINST  LOG ( PEP  CENT  RESIDUE  WEIGHT) 

c 

CALL  GP(PIB,WF,L.LS,M,JN,LW,LN,  A,  PLOT) 

ise 

WRITE  (6,3000) 

C _ _ 


c 

c 

LOOK  F00  FURTHER-  SETS  HF  DATA 

185 

3C0 

READ  (5 
I F( MORF 

,1300)  MORE 
."0.3)  GO  TO  1 

19C 

PO 

1000 

STOP 

FORMAT 

( 2X, A3, ?x, 8A6) 

1100 

FORMAT 

( 10X, A3,?X,8A6) 

1200 

FORMAT 

(  11, IX,  A4,  1 4, El  3.5, F 6.1, FI  3.5, F 6.1,  FI  3.5, F6. 7  ) 

1300 

FORMAT 

(ID 

1400 

FORMAT 

( 10X,  I3,4X,-3PFT.3, 5X.F6.3, 5X , 0PF6 . 3 , 2 ( 5 X , F6 . 3 ) ,2(5X,F6.4) ) 

1425 

FORMAT 

( / / 10X , 20H  AVERAGE  ACTIVATION  ENERGY  =  .-3PF6.3) 

143  5 

FORMAT 

( 1 C X , 1 7 H  AVERAGE  LOG  PREX,10X,2H=  ,F6.3) 

1 44 0  FORMAT! 10X, 34H30TH  FOR  10- 80  PERCENT  WEIGHT  LOSS) 

1  SCO  F  OR  MAT _ (  /J?  Xj_TH  WT  LF’SS ,  ?X  ,  SHE  A  (  KCAL  )  ,  3X  ,  9HST.  OE  VN.  ,  3X  ,  8HL0G  PPEX, 

.3X, 9HST.0EVN., ’X, 1CHAV.LOG  AFW , 2X , RHST.DE VN. , 2X, 1 1HL0G  RES.WT.) 

1700  F QRMAT  ( IPX, 1 8HL0G  R  ATr  VS  1/TEMP /IPX, 14HWE IGHT  LOSS  =  ,14) _ 

1 80 n  FORMAT  (  / 1  OX ,  1 1HPUN  10  NOS  ,9(A4,2H,  )) 

10 0  0  FORMAT  ( 1CX, 13HEPBQR  FOR  w  =,14,7HRUN  NO  ,A3,6H  PFAD  ,13, _ 

•OH  INST fad.) 

2DC 6  FORMAT  ( 1CX , 2  5  HL  FS  S  THAN  3  HEATING  PATES/1H1) _ 

30F0  FORMAT  (lHl') 

31 0  0  _F_OR  MAT  (1C  X  ,  3  2HACT  I  VAT  ION  ENERGY  VS  WEIGHT  LOSS) _ 

3?C0  FORMAT  (  1C-X,  22HPRF-FXP  VS  WEIGHT  LOSS) 

3  3  C  0  F  0 IRMA.  T  (  IPX  ,  46HAVFP  LOG  AF(w)  VS  LOG  PERCENT  WEIGHT  REMAINING) 

END 
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RATES  OF  WEIGHT  LOSS  FOR  A 
TYPICAL  TEFLON  EXPERIMENT 
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1CTF  19/03/68  TEFLON  (NEW  BAL.) 

CHART  SPEEC  (INCHES  PER  MINUTE )  -1.0000 
CHART  SPEEC  (FURLONGS  PER  FORTNIGHT!  -  2.5455 
OFFSET  (INCHES)  *  0.0625 
XRANGE  (INCHES)  -  36.60 

INITIAL  TEMP  (C)  -  418 

FINAL  TEMP  (Cl  -  625 

NC  OF  PTS  IN  CURVE  'IT  -  10 
TOTAL  NO  OF  POINTS  *  200 

MAX  TEMP  ERROR  -  2.616966 

TEMP  ROOT  MEAN  SQUARE  ERRCR  -  0.635928 

TEMP  POLY  COEFF 
0.000003 
-0.000230 
0.007C66 
-0.0872! 7 
5.973110 
418.330875 


HEIGHT  LOSS 

CWOT(NM) 

TPOLY 

TOER 

RTEMP 

1 

C.43704E— OC 

507.899 

0.57C00E 

01 

0. 12803E-02 

2 

0.76467E 

OC 

518.102 

0.57125E 

01 

0. 12638E-02 

3 

0 . 1053  IE 

01 

524.516 

0.57103E 

01 

0. 12536E-02 

4 

0.13932E 

01 

529.237 

0.57C40E 

01 

0.12463E-02 

5 

C. 17160E 

01 

532.890 

0.56964E 

01 

0. 12406E-02 

6 

C.19640E 

01 

535.986 

0.56884E 

01 

0.12359E— 02 

7 

C.22448E 

01 

538.693 

0.56803E 

01 

0. 12318E-02 

8 

C.25599E 

01 

54 1  uC67 

0.56724E 

01 

0. 12282E-02 

9 

0.29225E 

01 

543.174 

0.56649E 

01 

0.12250E-02 

10 

0.32235E 

01 

544.993 

0.56581E 

01 

0. 12223E-02 

11 

0.34049E 

01 

546.670 

0.56516E 

01 

0. 12198E-02 

12 

C.35286E 

01 

548.339 

0.56449E 

01 

0.12173E-02 

13 

C.37758E 

01 

549.891 

0. 56385E 

01 

0. 12150E-02 

14 

C.40299E 

01 

551.333 

0.56326E 

01 

0. 12129E-02 

15 

C.42162E 

01 

552.693 

0.56269E 

01 

0. 12109E-02 

16 

0.44440E 

01 

553.998 

0.56214E 

01 

0. 12090E— 02 

17 

0.46574E 

01 

555.237 

0.56162E 

01 

0. 12072E-02 

18 

0.500UE 

01 

556.430 

0.56113E 

01 

0. 12054E-02 

19 

C.53156E 

01 

557.516 

0.56C68E 

01 

0. 12038E-02 

20 

0.55268E 

01 

558.527 

0.56C27E 

01 

0. 12024E-02 

21 

C.56799E 

01 

559.508 

0.55988E 

01 

0. 12010E-02 

22 

C.57782E 

01 

560.473 

0.55950E 

01 

0. 11996E-02 

23 

0.58077E 

01 

561.465 

0.55912E 

01 

0. 1 1901E-02 

24 

0.59939E 

01 

562.406 

0.55877E 

01 

0. 1 1968E-02 

25 

C.61520E 

01 

563.326 

0.55843E 

01 

0. 1 1955E-02 

26 

C.63L33E 

01 

564.223 

0.55812E 

01 

0. 1 1942E-02 

27 

0.65319E 

01 

565.113 

0. 5 5782 E 

01 

0. 1 1929E-02 

28 

C.68436E 

01 

565.958 

0.55754E 

01 

0. 11917E-02 

29 

C.70771E 

01 

566.750 

0.55729E 

01 

0. 1 1906E-02 

30 

0.73456E 

01 

567.515 

0.55707c 

01 

0. 1 1895E-02 

31 

0.74898E 

01 

568.251 

0.55686E 

01 

0. 1 1885E-02 

32 

0.75789E 

01 

568.978 

0.55666E 

01 

0. 1 1875E-02 

33 

0.76239E 

01 

569.701 

0.55648E 

01 

0. 11864E-02 

34 

C.76522E 

01 

570.446 

0.55631E 

01 

0. 1 1854E-02 

35 

C.77316E 

01 

571.172 

0.55615E 

01 

0. 1 IB44E-02 

36 

0.78633E 

01 

571.894 

0.55601E 

01 

0. 11834E-02 

37 

0 . 80040E 

01 

572.598 

0.55588E 

01 

0. 1 l 824E-02 

38 

0.82437E 

01 

573.284 

0.55577E 

01 

0. 1 1814E-0? 

39 

0.83949E 

01 

573.950 

0.55568E 

01 

0.11805E-0? 
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40 

0.85006E  01 

574.597 

0.55560E  01 

0. 11796E-02 

41 

0.85905E  01 

575.245 

0.55554E  01 

0. 11787E-02 

42 

C.87186E  01 

575.883 

0.55549E  01 

0.U778E-02 

43 

0.87468E  01 

576.515 

0.55545E  01 

0. 11769E— 02 

44 

0.87688E  01 

577.151 

0.55543E  01 

0.  11760E-02 

45 

0.88183E  01 

577.783 

0.55543E  01 

0.11752E-02 

46 

0.88743E  01 

578.419 

0.55543E  01 

0.  11743E-02 

47 

0.91338E  01 

579.047 

0.55546E  01 

0. 11734E-02 

43 

0.92862E  01 

579.648 

0.55550E  01 

0.11726E-02 

49 

C.94210E  01 

580.236 

0.55555E  01 

0.U718E-02 

50 

C.95513E  01 

580.817 

0.55562E  01 

0.11710E-02 

51 

0.97158E  01 

581.385 

0.55570E  01 

0.11702E-02 

52 

C.97623E  01 

581.946 

0.55579E  01 

0.1 1694E-02 

53 

0.97623E  01 

582.510 

0.55590E  01 

0. U687E-02 

54 

0.97445E  01 

583.086 

0.55602E  01 

0. 116796-02 

55 

0.97178E  01 

583.674 

0.55617E  01 

0.116716-02 

56 

C.99691E  01 

584.248 

0.55632E  01 

0.11663E-02 

57 

C.10038E  02 

584.802 

0.556495  01 

0. 11656E-02 

58 

0.10131E  02 

585.351 

0.55668E  01 

0.11643E-02 

59 

C.10208E  02 

585.896 

0.55687E  01 

0. 116416-02 

60 

C.10281E  02 

586.432 

0.55708E  01 

0.  U633E-02 

61 

C.10400E  02 

586.964 

0.55731E  01 

0. 11626E-02 

62 

C.10414E  02 

587.493 

0.55755E  01 

0.11619E-02 

63 

C.10376E  02 

588.031 

0.5578CE  01 

0. 11612E-02 

64 

C.10310E  02 

588.572 

0.55808E  Cl 

0.11605E-02 

65 

0.10204E  02 

589.117 

0.55838E  01 

0. 11597E-02 

66 

C.10352E  02 

589.667 

0.55870E  01 

0. 11590E-02 

67 

0.10344E  02 

590.210 

0.55903E  01 

0. 11583E-02 

68 

0.10358E  02 

590.751 

0.55938E  01 

0. 1 1575E-02 

69 

0.10391E  02 

591.290 

0.55975E  01 

0.  H568E-02 

70 

0.10387E  02 

591.825 

0.56C13E  01 

0. 11561E-02 

71 

C.10423E  02 

592.364 

0.56054E  01 

0. 11554E-02 

72 

0.10433E  02 

592.894 

0.56096E  01 

0. 11547E-02 

73 

0. 10321E  02 

593.423 

0.56140E  Cl 

0. 115406-02 

74 

0.10237E  02 

593.962 

0.5&186E  01 

0.11532E-02 

75 

C.10099E  02 

594.509 

0.56235E  01 

0. 1 1525E-02 

76 

0. 100 11E  02 

595.089 

0.56290E  01 

0. 1 1517E-02 

77 

0.10151E  02 

595.659 

0.56346E  01 

0. 1 1510E-02 

78 

0.10124E  02 

596.219 

0.56403E  01 

0. 1 1502E— 02 

79 

0.10033E  02 

596.774 

0.56463E  01 

0. 1 1495E-02 

80 

C.99414E  01 

597.328 

0.56524E  01 

0. 1 1488E-02 

81 

0.97937E  01 

597.889 

0.56589E  01 

0.  11480E-02 

82 

0.96724E  01 

598.456 

0.56657E  01 

0.  11473E-02 

83 

0.94010E  01 

599.041 

0.567295  01 

0.  1 1465E-02 

84 

C.90958E  01 

599.649 

0.56808E  01 

0.11457E-02 

85 

0.86745E  01 

600.311 

0.56897E  01 

0.  U449E-02 

86 

0.84766E  01 

600.964 

0.5698SE  01 

0.11440E-02 

87 

0.82320E  01 

601.653 

0.57C88E  01 

0.11431E-02 

88 

0.79915E  01 

602.368 

0.57196E  01 

0.  1 1422E-02 

89 

0.77573E  01 

603.100 

0.573UE  01 

0.11412E-02 

90 

0.75096E  01 

603.846 

0.57433E  01 

0.1 14026-02 

91 

0. 72152E  01 

604.618 

0.57565E  01 

0.  1 1392E-02 

92 

C.68785E  01 

605.430 

0.57709E  01 

0.11382E-02 

93 

0.64873E  01 

606.282 

0.57866E  01 

0.  1 1371E-02 

94 

C.59624E  01 

607.219 

0.68C47E  01 

0.  1 1359E-02 

95 

0.54449E  01 

608.232 

O.50252E  01 

0.  11346E-02 

96 

C.48352E  01 

609.365 

0.58493E  01 

0. 1 133  IE-02 

97 

0.41379E  01 

610.678 

0.58787E  01 

0. 1 1314E-02 

98 

0.33280E  01 

612.262 

0.59164E  01 

0. 1 1294E-02 

99 

0.23131E  01 

614.390 

0.597096  01 

0.11267E-02 

AVERAGE  TEMP  DERIVATIVE  • 

0.56306E 

01 
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teflon  degradation 


H 

.A.  *  HEATING  RATE 

IN  DEGREES  C 

PER  HOUR 

RATE  CF  WT. 

LOSS  IN 

PERCENT  PER 

m?nute 

—  . - 

-  . 

I .Q./H.R. 

4TF/ 

85 

9TF/145 

11TF/44 

1 2  TF  / 

55 

PERCENT 

RATE  OF 

TEMP. 

RATE  OF 

TEMP. 

RATE  OF 

TEMP. 

RATE  OF 

TEMP. 

NT. LOSS 

NT. LOSS 

OEG  C 

NT  .LOSS 

OEG  C 

NT. LOSS 

OEG  C 

MT .LOSS 

OEG  C 

I 

0.074 

486.5 

0.207 

494.9 

C.069 

480.8 

0.078 

479.9 

2 

0.210 

497.6 

0.390 

505.3 

0.138 

489.4 

0.154 

489.4 

3 

0.312 

503.4 

0.553 

511.2 

0.194 

494.5 

0.212 

494.6 

4 

0.418 

507.6 

0.708 

515.6 

0.243 

498.3 

0.265 

498.5 

5 

0.526 

510.8 

0.857 

519.0 

C.307 

501.3 

0.323 

501.7 

6 

0.595 

513.4 

0.998 

521.9 

0.349 

503.6 

0.383 

5G4.3 

7 

0.692 

515.8 

1.134 

524.4 

C  .403 

505.8 

0.433 

506.6 

8 

0.776 

517.8 

1.286 

526.6 

0.451 

507.6 

0.504 

508.7 

9 

0.825 

519.6 

1.412 

528.5 

C.509 

509.2 

0.594 

510.4 

10 

0.910 

521.3 

1.590 

530.2 

0.563 

510.7 

0.634 

511.9 

u 

0.987 

522.8 

1.722 

531.8 

0.598 

512.0 

0.637 

513.4 

12 

1.091 

524.3 

1.811 

533.2 

C.640 

513.3 

0.688 

514.8 

13 

1.135 

525.5 

1.932 

534.6 

C.688 

514.4 

0.732 

516.1 

14 

1.239 

526.7 

2.054 

535.8 

0.720 

515.5 

0.789 

517.3 

15.  _ 

1.285 

527.8 

2.193 

537.0 

0.773 

516.5 

0.838 

518.4 

16 

1.363 

529.0 

2.327 

538.1 

0.805 

517.5 

0.882 

519.5 

IT. 

1.451 

530.0 

2.398 

539.2 

0.840 

518.4 

0.921 

520.6 

18 

1.532 

530.9 

2.513 

540.3 

0.868 

519.3 

0.981 

521.6 

.19 

1.6  22 

531.9 

2.678  . 

541.2 

0.948 

52C.1 

1.028 

522.5 

20 

1.704 

532.7 

2.813 

542.1 

C.989 

520.9 

1.066 

523.4 

21 

1.738 

533.6 

2.873 

543.0 

1.009 

521.6 

1.107 

524.2 

22 

1.793 

534.4 

2.878 

543.9 

1.023 

522.3 

1.146 

525.1 

2? 

1.887 

535.2 

2.997 

544.7 

1.021 

523.1 

1.207 

525.9 

24 

1.988 

535.9 

3.119 

545.5 

1.054 

523.6 

1.251 

526.6 

25 

2.053 

536.6 

3.155 

546.3 

1.086 

524.5 

1.301 

527.4 

26 

2.122 

537.3 

3.250 

547.0 

1.112 

525.2 

1.311 

528.1 

2T 

2.176 

538.0 

3.359 

547.0 

1.145 

525.8 

1.351 

528. e 

28 

2.207 

538.6 

3.455 

548.5 

1.159 

526.5 

1.415 

529.5 

29 

2.236 

539.2 

3.674 

549,2 

1.193 

527.1 

1.443 

530.1 

30 

2.255 

539.9 

3.812 

549. e 

1.217 

527.7 

1.468 

530.7 

31 

2.277 

540.5 

3.856 

550.5 

1.246 

528.3 

1.484 

531.4 

32 

2.327 

541.1 

3.726 

551.1* 

1.282 

528.9 

1.526 

532.0 

JJ 

2.384 

541.7 

3.888 

551.7 

1.313 

529.4 

1.578 

532.6 

34 

2.459 

542.3 

3.924 

552.3 

1.334 

530.0 

1.618 

533.2 

35 

j2*537 

342.9 

4.009 

552.9 

1.367 

530.5 

1.670 

533.7 

36 

2.608 

543.5 

4.060 

553.5 

1.393 

531.1 

1.672 

534.3 

U 

2.688 

544. Q 

4.108  _ 

514.1 

1.413 

531.6 

1.710 

534.9 

38 

2.754 

544.5 

4.222 

554.7 

1.447 

532.1 

1. 749 

535.4 

39 

£.026 

3*5,0 

4.206 

555.2 

L.469 

532.6 

1.796 

535.9 

40 

2.875 

545.5 

4.292 

555.8 

1.489 

533.1 

1.788 

536.4 

51 

1,912 

546.0 

4.366 

556.4 

1.504 

533.6 

1.792 

536.9 

42 

2.940 

546.5 

4.400 

556.9 

1.525 

534.0 

1.827 

537.5 

.  .43. 

2,989 

546.9 

4.493 

557.4 

1.538 

534.5 

1.845 

538.0 

44 

3.014 

547.4 

4.537 

558.0 

1.542 

535.0 

1.846 

538.5 

_.45 

3.-025  -  _ 

547,9 

4.603 

558.5 

1.554 

535.4 

1.854 

539.0 

46 

3.032 

548.3 

4.628 

559.0 

1.564 

535.9 

1.827 

539.5 

41 

3.035 

548.8 

4.602 

559.5 

1.564 

536.4 

1.861 

540.0 

48 

3.048 

549.3 

4.706 

560.0 

1.564 

536.8 

1.876 

540.5 

49 

3013 

549.7 

4.783 

560.5 

1.564 

537.3 

1.865 

541.0 

50 

3.114 

550.2 

4.896 

561.0 

1.549 

537.7 

1.856 

541.5 
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PtuCeNT 

RAT  £  OF 

TEMP. 

RATE  OF 

TEMP. 

RATE  OF 

TEMP. 

RATE  OF 

TEMP. 

WT.LOSS 

NT. LOSS 

DEG  C 

MT  .LOSS 

DfG  C 

WT.LOSS 

&fITC 

If  .1353 

Deg  c 

51 

3.128 

550.6 

4.917 

561.5 

1.542 

538.2 

1.869 

542.0 

52 

3.106 

551.1 

47903 

561.9 

i:535 

1.874 

54274~ 

53 

3.108 

551.5 

4.951 

562.4 

1.521 

539.2 

1.858 

542.9 

54 

3.106 

552.0 

4.978 

562.9 

1.5T3 

339.6 

1.838 

34374“ 

55 

3.125 

552.4 

5.055 

563.4 

1.517 

540.1 

1.832 

543.9 

56 

3.139 

552.9 

5.041 

563.8 

1.513 

54576 

1.827 

33474 

57 

3.146 

553.3 

5.047 

564.3 

1.514 

541. C 

1.898 

545.0 

58 

3.111 

553.8 

5.120 

564.8""" 

~ 1.506  " 

541.3 

1.925 

543.4 

59 

3.116 

554.2 

5.144 

565.2 

1.504 

542.0 

1.930 

545.9 

60 

3.100 

554  c  7 

5.196 

565.7 

1.495 

542.5 

1.890 

546.4 

61 

3.136 

555.1 

5.185 

566.1 

1.480 

543.0 

1.810 

546.9 

62 

3.125 

555.6 

5.150 

566.6 

1.463 

543.  f 

1.816 

547.4 

63 

3.113 

556.0 

5.190 

567.0 

1.448 

543.9 

1.776 

547.9 

64 

3.086 

556.4 

5.179 

56  775 

1.434T 

344.4 

T7752 

34874 

65 

3.052 

556.9 

5.215 

567.9 

1.425 

544.9 

1.732 

549.0 

66 

3.014 

557.4 

5.207 

568.4 

1.429 

545.4 

1.763 

34973~~ 

67 

2.980 

557.8 

5.138 

568.3 

1.440 

546.0 

1.779 

550.0 

68 

2.948 

558.3 

5.221 

569.3 

1.437 

546.4 

1.767 

550.5 

69 

2.935 

558.8 

5.245 

569.7 

1.419 

546.9 

1.714 

551.0 

70 

"2.920 

559.3 

5.2/0 

570.1 

1.390 

54T.4 

i 

"35T76 

71 

2.895 

559.8 

5.268 

570.6 

1-379 

547.9 

1.707 

*52.1 

7? 

2.862 

560.2 

5.175 

571.0 

1.357 

548.5 

1.716 

552.6 

73 

2.890 

560.7 

5.108 

571.5 

1.348 

549.0 

1.721 

553.2 

74 

2.880 

561.2 

5.065 

571.9 

1.361 

549.5 

1.696 

553.7 

75 

2.861 

561.7 

4.968 

572.4 

1.368 

550.1 

1.686 

554.2 

76 

2.846 

562,2 

4.933 

572.9 

1.368 

550.6 

1.683 

554.8 

77 

2.829 

562.7 

4.862 

573.4 

1.368 

551.1 

1.688 

555.3 

78 

2.779 

563.2 

4.914 

573.8 

1.368 

551.6 

1.743 

555.0 

79 

2.721 

563.7 

4.949 

574.3 

1.368 

552.1 

1.761 

556.4 

80 

2.686 

564.2 

4.931 

574.8 

1.368 

552.7 

1.739 

356.9 

61 

2.631 

564.7 

4.830 

575.2 

1.368 

553.2 

1.730 

557.4 

82 

2.577 

565.3 

4.666 

575.7 

1.368 

553.7 

1.680 

55773 

e3 

2.508 

565.8 

4.476 

576.3 

1.368 

554.2 

1.668 

558.5 

84 

2.514 

566.4 

4.350 

576.8 

1.368 

554.7 

1.637 

559.0 

85 

2.475 

567.0 

4.204 

577.3 

1.354 

555.2 

1.627 

559.6 

86 

2.443 

56  7.5 

4.093 

T?7.~<r 

1.354 

555.8 

1.632 

560.1 

87 

2.411 

568.1 

3.996 

578.5 

1.347 

556.3 

1.614 

560.7 

88 

2.349 

568.7 

3.851 

579.1 

17326 

556.8 

“17378 

36172 

89 

2.278 

569.3 

3.734 

579.7 

1.303 

557.4 

1.539 

561.8 

SO 

2.218 

569.9 

3.574 

580.4 

1.272 

357.9 

17489 

36274 

91 

2.145 

570.5 

3.442 

581.0 

1.235 

558.5 

1.420 

563.0 

92 

2.049 

571.3 

3.188 

581.8 

1.196 

559.1 

1.346 

363.7 

93 

2.009 

571.9 

3.016 

582.5 

1.158 

559.7 

1.255 

564.3 

94 

1.922 

572.6 

2.845 

583.4 

1.104 

5603“ 

17168 

56371" 

95 

1.818 

573.4 

2.652 

584.2 

1.028 

560.9 

1.059 

565.9 

96 

1.664 

574.1 

2.459 

585.2 

0.923 

561.6 

0.946 

566.8 

97 

1.454 

575.1 

2.180 

586.2 

0.794 

562.5 

0.817 

567.8 

98 

1.181 

576.1 

1.783 

587.5 

0.608 

563.5 

0.649 

569.0 

99 

0.827 

577.6 

1.115 

589.2 

0.392 

564.9 

0.429 

570.6 
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XXX  FRIEDMANS  * FT HDD  FQ8  EA«-A 


RUN  10  NOS 


WT  LOSS  EAIKCALI 


69.464 

'Tf7355 

67.468 

61.834 

&T.T75 

66.386 


65.037 
TT4.998 
64.214 
~SV.  576 
64.450 
6474J6 
67.620 
64.320 
63.404 
62.178 
63.481 
64.328 
63.056 
64.061 
^7.845 
63.  ■j7't 
65.678 
66. 5  <»5 
66.057 
62.574 
63.252 
6?,930" 
62.409 
62.726 


9TF,  1ITF,  12TF, 


ST.OEVN. 

— tttbt 

0.758 

- 1.496 

0.842 

~"T77ET 

2.756 


1.741 

“177576“ 
0.724 
’  37998' ~ 
1.656 
"T75T5 
0.962 
"  C.TTT ' 

1.1  84 

r.?ro“ 

1.137 
0.61 1 
0.486 
0.6  85" 

0.95’. 

0.951 ~ 
0.807 
1.046 
0.670  _ 
0.8 27' 
C.  ’’04 
0.384  " 
0.564 
" 0.685 
C.806 
0.738 


LOG  PPEX 
19.591' 
19.092 
"TW 
18.536 
”"17.767“ 
16.991 

"  177479 . 

18.245 

17.877 

17.864 

17.650 

177374 

17.710 

T77W~ 

17,492 

1 7.68  6”' 

17.427 

17.085 

17.546 

17.669 

17.321 

17.588 

1 7.527 

17. 55* 

1 8. C19 
“18.261 
1 8.1 08 
1  77161" 

17.345 
If.  25  5 
17,115 
1T7T9  4 


ST.OEVN.  A V. LOG  AFW  ST.OEVN.  LOG  RES.WT. 


C.453 


0.501 


1.631 

1.488 

1.026 

“TTFtrr 

C.426 

075T7" 

C.912 

0.900 

0.562 

07T56” 

C  .690 

0.704 

C.662 

C.355 

0.283 

C.307 

C  .6  51 

0.55i 

C  .467 

C  .604" 

_ C  .  3  87 

0.477 
C  .406 
0.221 
C  .325 
"C  .394 
C.463 
C74T9 


19.058 

19.057 

19.068 

19.071 

19.060 

"T97T55TT' 

19.C59 

"17.35?" 
19.054 
19.051“ 
19.046 
T57045~ 
19.051 
T7.T50' 
19.042 
19.031  ' 
15.028 
19.028 
19.023 

19  ToTT 

19.014 

19.009 

19.008 

19.002" 

!0.993 

18.984 

13.985 
18. 980 

18.979 

"TT797? 


0.0065 


1.5777 


0.0096  1.9685 


0.0351  1.9590 


0,01 93 
D73T5F — 
0.0156 
3. 01 67 
0.0223 

■"373734 

0.0184 

""373173 

0.0219 

"073705“- 
0.0224 
0.924 j 
0.0183 
0.0177 
0.0226 
"370195 
0.0198 
575204" 
0.0137 
0.0"I73“ 
0.0127 

'  c  .3578"' 
7U021C 
0.0222" 
0.0244 

0.3H3 


1.9494 


1 .  9  395 


1.9294 

“173743“ 

1.9191 

“175173“ 
1.9C85 
“1.9351“ 
1.9576 
1.3521 
1.8865 
1.8 8C8 
1.8751 
“TTTfO? 
1.8633 
1775T3 
1,8£1? 

"  1.5451 
1.8  388 
'  173325“ 
1.8  261 
"  175198 
1.°  129 
1".  .67 


37 

63.605 

0.300 

17.43C 

C.459 

18.97C 

C  .0207 

1.7593 

39 

67. 6" 3 

1.059 

17.159 

CTOOfT 

13.566 

T73T47 

1.7=24  '"' 

39 

62.500 

1.287 

17.1 3C 

C.733 

18.964 

0.0262 

1 . 7  653 

40 

62.142 

1.485 

”  17.026 

r73?r 

IT.  955 

3.3281 

1.7782 

41 

62.613 

1.515 

17.146 

f  .866 

18.947 

0.0271 

1.7709 

42 

62.371 

1.4  35 

1  6  799  b 

~‘f  7812 

“  13.94? 

0.3 fa 5 " 

1.7e34 

43 

62.93? 

1.553 

17.223 

C.989 

19.937 

0.0265 

1.7559 

44 

67.25? 

17655 

17.TCT 

“C.?46 

18T " 

3.0263 

1.7482  " 

45 

63.56? 

1.490 

17.376 

C  .  8 56 

18.920 

0.0244 

1.7404 

46 

64.11? 

1.677 

17.513'  — 

" C.95? 

13.9(57 

3.0242  “ 

"T.7  324 - 

47 

63.415 

1.607 

17.115 

C.C17 

18.899 

0.0255 

1.7243 

48 

64. ?c5 

l.?78 

3 

*7 

1 

l 

i 

C  v  786 

19.392 

'  5.0217 

1.7160 

49 

65.801 

1.797 

17.941 

1.024 

18.885 

0.0217 

1.7C76 

50 

67.473 

1.644 

"  '  1 A7T7? 

C.  93f 

- 18.875 

373175 

1 .0  530 

51 

67. ??4 

1.  741 

i8.4ee 

C.-9I 

10.8t6 

0.0183 

1.6502 

52 

63.027 

1.583 

18.5C5 

0.901 

18.854 

0.0164 

1.6612 

53 

69.371 

!  .665 

18.94C 

0.947 

13.943 

0.016° 

1.6721 

54 

69. 9?l 

1.635 

18.9PP 

C.958 

13.831 

0.0168 

1.6629 

55 

71.029 

1.769 

1°.271 

C.9°9 

18.323 

0.0186 

1.6632 

56 

71.454 

1.906 

19.372 

l  .3C2 

13.512 

o.ofol 

1.6435 

57 

70.581 

l.«19 

19.134 

1.032 

18.806 

0.0185 

1.6  335 

58 

71.095 

1.513 

19.261 

C.858 

18. 796 

0.0159 

l .6232 

59 

71.811 

1.513 

19.441 

0.369 

18.787 

0. Jl 73 
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60 

73.201 

1.384 

19.795 

0.784 

18.774 

0.0135 

1.6C21 

61 

75.168 

1.968 

20.301 

1.114 

18.759 

0.0272 

1.5511 

62 

74.903 

1.965 

20.241 

1.112 

18.746 

0.0267 

1.5798 

63 

76.869 

2.105 

20.723 

1.190 

18.733 

3.0319 

1.5682 

64 

77.521 

2.216 

20.880 

1.252 

18.719 

0.9342 

1.5  563 

65 

78.950 

2.124 

21.242 

1 .200 

18.705 

0.0371 

1.5441 

66 

78.260 

1.632 

21.048 

0.921 

19.695 

0.0327 

1.5215 

67 

77.474 

1.528 

20.828 

r.862 

13.684 

0.0296 

1.5185 

68 

78.446 

1.266 

21.072 

C  .  7 1 4 

18.673 

0.0313 

1.5C51 

69 

80.437 

1.562 

21.68C 

0.88Q 

18.658 

0.0381 

1.4514 

70 

82.296 

1.469 

22.054 

C  .828 

1 8 .644 

0.0431 

1.4771 

71 

82.603 

1.260 

22.120 

C.709 

18.631 

0.3433 

1.462*. 

72 

82.802 

1.074 

22.155 

0.604 

18.616 

0.0431 

1.4472 

73 

82.576 

1.515 

22.081 

C.851 

18.604 

0.0434 

1.4  214 

74 

82.425 

1.599 

22.028 

0.898 

18.592 

0.0430 

1.4150 

75 

81.455 

1.733 

21.756 

0.973 

18.577 

0.0405 

1.357? 

76 

81.238 

1,328 

21.685 

1.026 

18.564 

9.0401 

1.3802 

77 

80.322 

1.756 

21.430 

e.984 

18.552 

0.0373 

1.3617 

78 

80.557 

0.855 

21.481 

c.47C 

1«.543 

0.0350 

1.3424 

79 

80.596 

0.447 

21.478 

C.2  5C 

18.532 

0.0342 

1.3222 

80 

80.958 

0. 

21.557 

0. 

18.513 

0.0351 

1 »  3  CIO 

81 

80.003 

0.390 

21.290 

C.218 

19.502 

0.0322 

1.2788 

82 

78.423 

0.804 

20.355 

C  .450 

18.482 

0.0282 

1.2553 

83 

75.586 

0.929 

20.0=2 

0.519 

>C 

* 

• 

<r 

J  .  0234 

1.2  204 

84 

74.344 

1.073 

19.751 

9  •  5  c  9 

13.445 

0.0179 

1.2C41 

85 

72.872 

0.996 

19.347 

C.556 

18.426 

3.0139 

1 . 1 731 

86 

70.995 

0.730 

18.841 

C  .407 

18.410 

0.0093 

1.1431 

87 

69.734 
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